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INTRODUCTION 

Research  into  the  behavior  of  projectiles  pene- 
trating into  seafloor  materials  was  completed  recently 
by  the  Civil  Engineering  Laboratory  (CEL),  Naval 
Construction  Battalion  Center,  Port  lluencme, 
California.  The  results  of  this  study  are  directly 
related  to  the  design  of  direct-embedment  anchors. 
Embedment  anchor  systems  have  been  developed 
(Smith  1971,  and  Taylor  and  Beard,  1973)  with  rea- 
sonable success  by  using  existing  relationships  and 
exercising  a good  deal  of  judgment.  However,  it  was 
recognized  that  substantial  time  and  dollar  savings 
could  be  realized  in  future  programs  if  designers  had 
available  a method  for  making  more  accurate  and 
reliable  predictions  of  embedment  depth  for  a can- 
didate embedment  anchor  system.  This  report  deals 
with  the  development  of  such  a method. 

In  addition  to  the  design  of  direct  embedment 
anchors,  the  behavior  of  projectiles  penetrating  into 
soils  and  cemented  soil  materials  is  an  integral  part  of 
several  other  types  of  engineering  endeavor.  Piles 
resist  abrupt  settlement  according  to  the  same  laws  of 
mechanics  governing  penetration.  Although  the 
technology  of  pile  design  is  well-developed,  it  is 
believed  that  this  research  can  contribute  to  the 
development  of  more  accurate  pile  design  techniques 
for  deep  piles  in  soft  seafloor  soils.  Also,  the  applica- 
tion of  penetration  behavior  in  the  design  of  dynamic 
penetrometers  and  in  the  utilization  of  dynamic 
penetration  data  to  evaluate  the  engineering  proper- 
ties of  soils  is  important.  Work  presently  underway  at 
CEL  to  develop  an  optimum  expedient  penetrometer 
for  deep  seafloor  soils  is  fully  utilizing  the  results 
reported  herein. 

A further  application  of  penetration  technology 
is  in  the  designing  of  better  coring  equipment.  This 
application  is  currently  under  investigation  at  CEL. 

A final  category  of  seafloor  engineering 
benefiting  from  improved  penetration  technology  is 
the  estimation  of  the  depth  of  embedment  of 
salvageable  objects.  The  estimated  portion  of  an 


object  visible  above  the  seafloor  is  important  in  plan- 
ning a search,  and  the  embedment  depth  should  be 
considered  in  preselecting  gear  for  grappling, 
grabbing,  securing,  and  lifting.  For  most  salvage 
efforts,  embedment  is  relatively  shallow,  when 
compared  to  the  deep  penetration  of  primary  concern 
in  the  present  work.  However,  the  basic  conclusions 
of  the  present  work  relating  to  the  applicability  of 
conventional  bearing  pressure  theory  and  dynamic 
soil  strength  can  be  generalized  readily  to  the  shallow 
penetration  case. 

BACKGROUND 

The  problem  of  penetration  has  been  considered 
analytically  in  the  development  of  propellant- 
launched,  penetrating  anchors.  Christians  and 
Mcisburger  (1967)  relied  upon  an  empirically  derived, 
specific  form  of  the  Poncelct  Equation.  Others, 
including  Aerojet  (1968),  have  utilized  simpler 
expressions  involving  a constant  resistance  to  pene- 
tration. 

An  exhaustive  assessment  of  the  state  of  the  art 
in  penetration  prediction  was  performed  by  Schmid 
(1969).  He  recommended  integrating  equations  of 
motion  incorporating  basic  material  laws  (static  and 
viscous  resistance,  for  example)  in  order  to  obtain 
relationships  for  predicting  the  penetration  behavior 
of  free-falling  objects  into  soft  seafloor  materials. 
However,  the  choice  of  which  material  resistance  law 
was  left  to  the  user. 

An  extensive  study  of  high-speed  penetration 
into  earth  materials  was  performed  by  Sandia  Corpor- 
ation; the  results  were  summarized  by  Young  (1969). 
Most  of  the  Sandia  work  was  done  with  rocks  and 
harder  soils,  although  a few  tests  were  conducted  in 
soft  muds.  Young  recommended  the  use  of  an  empiri- 
cal nomograph  together  with  empirically  obtained 
soil  parameters  not  easily  related  to  conventional  soil 
test  data. 
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Field  tests  were  conducted  by  Migliore  and  Lee 
(1971)  at  CEL  on  an  instrumented  penetrometer/ 
corer  in  a terrigenous  mud  in  1 30  feet  of  water  in  the 
Santa  Barbara  Channel  to  acquire  data  for  use  in 
evaluating  proposed  penetration  prediction  tech- 
niques. Data  from  these  tests  showed  substantial 
variation  from  available  predictive  techniques.  A 
regression-analysis-derived  procedure  was  recom- 
mended for  making  predictions  in  soils  similar  to  the 
test  target,  but  generalization  to  other  soils  based 
upon  soil  test  data  or  other  soil  property  information 
was  left  in  doubt. 

A recurring  problem  in  the  works  cited,  as  well 
as  in  other  works  (see  the  references  given  kv  Schmid, 
1969,  for  example),  is  the  difficulty  encountered  in 
relating  theories  to  actual  soil  properties  except 
empirically  for  each  new  soil  of  concern.  A better 
understanding  was  needed  of  how  the  physical  and 
mechanical  properties  of  soils  as  determined  from 
conventional  soil  mechanics  testing  procedures  could 
be  incorporated  with  physical  laws  into  a general  rela- 
tionship for  penctrator  performance. 

APPROACH 

The  approach  taken  in  this  work  involved  pri- 
marily the  use  of  scale-model  tests  to  determine 
empirically  the  relationships  between  soil  properties 
and  penctrator  performance.  Test  data  were  analyzed 
within  a framework  that  was  fairly  rigidly  predeter- 
mined by  the  laws  of  motion,  the  theories  of 
plasticity  and  fluid  flow,  and  conventional  practice  in 
the  areas  of  soil  testing  and  foundation  bearing 
capacity  assessment.  Data  from  field  tests  and  other 
sources  were  used  to  substantiate,  refine,  and  verify 
the  scaling  up  to  large  size  of  the  relationships  derived 
from  the  model  test  program.  User  recommendations 
were  based  upon  these  relationships  with  constants 
derived  from  the  testing,  from  accumulated  data  on 
seafloor  soil  properties,  and  from  conventional  fluid 
mechanics  experience, 

MODEL  TESTS 
Test  Procedure 

Each  test  consisted  of  firing  a penctrator  into  a 
target  while  taking  high-speed  movies  of  the  penc- 


trator and  target  and  monitoring  the  acceleration  of 
the  penctrator. 

Target  Materials.  Materials  used  in  the  model 
tests  were  (1)  a clayey  silt  from  Scat  Beach,  Cali- 
fornia, in  soft  and  medium-soft  consistencies,  (2)  the 
same  silt  together  with  poftland  cement  in  three 
different  formulations  of  soil-ccmcnt  having  different 
strengths,  and  (3)  port  lend  cement  concretes  made 
with  expanded  shale  aggregate  in  soft  and  hard  for- 
mulations and  with  a river  sand  aggregate  in  a formu- 
lation of  medium  strength. 

Silt  targets  were  constructed  by  placing  the  soil, 
which  was  in  a pasty  consistency,  into  a specially 
constructed  tank.  The  paste  was  prepared  in  a 
planetary  mixer  under  a vacuum  of  approximately  29 
inches  of  mercury  in  order  to  minimize  the  amount 
of  free  air  left  in  the  soil  after  placement.  It  was  felt 
that  this  procedure  would  produce  a uniform  target 
in  which  the  material  was  essentially  saturated, 
thereby  giving  a true  indication  of  the  variations  in 
soil  shearing  resistance  with  velocity.  However,  it  was 
impossible  to  remove  all  the  air  bubbles;  conse- 
quently, the  material  was  subject  to  some  local  com- 
pression rather  than  being  incompressible  as  had  been 
desired.  Although  the  small  air  bubbles  present 
undoubtedly  compressed  during  projectile  penetra- 
tions, it  is  felt  that  the  volume  accommodated  by 
compressing  bubbles  was  substantially  less  than  that 
displaced  by  the  advancing  projectile.  Therefore,  the 
observed  soil  deformation  patterns  and  penetration 
resistances  remained  indicative  of  the  behavior  that 
would  occur  in  a natural,  saturated  sediment. 

Each  prepared  target  was  used  repeatedly  for 
several  penetration  tests.  Three  half-section  tests  were 
conducted  at  different  locations  adjacent  to  the 
acrylic  plastic  face  of  the  test  tank.  Subsequently, 
several  full-section  tests  with  smaller  projectiles  were 
conducted  at  other  locations  in  the  remainder  of  the 
soil.  Care  was  taken  to  space  the  test  locations  several 
inches  apart  in  order  to  minimize  the  effects  on  a new 
test  of  a hole  left  from  a previous  test. 

Each  silt  target  was  prepared  by  placing  the  soil 
with  the  tank  in  a horizontal  orientation.  After  filling 
and  scrccding,  the  surface  to  be  covered  with  the 
acrylic  plastic  face  was  covered  with  a thin  layer  of 
graphite  and  contrasting  talcum-powder  lines  spaced 
on  a 2-inch-square  grid.  The  graphite  was  used  to 
reduce  frictional  contact  with  the  plastic  in  order  to 
give  an  accurate  picture  of  the  soil  deformations  in 
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Table  1.  Properties  of  Ccmcnted-Target  Materials 


Aggregate 

Strength 

Designation 

/Weight  Free\ 
\ Water  / 

Weight  Cement 
/Weight  Dry  \ 
\ Aggregate  / 

l Average  Shear 
Strength 
dcgf/cm2.)1* 

Weight  Cement 

Silt 

soft 

12.2 

0.06 

1.2 

Silt 

medium 

5.8 

0.13 

5.4 

Silt 

hard 

3.1 

0:26 

22 

Expanded  shale 

soft 

i.o 

0.26 

40 

Expanded  shale 

hard 

0.5 

0.58 

172 

Fine  sand 

- 

1.0 

0.31 

88 

a Specimen  and  target  aged  at  least  7 days;  strengths  measured  shortly  after 
penetration  testing,  softer  materials  by  vane  shear  and  harder  materials  by 
unconfined  and  confined  compression  tests. 


the  half-space.  Figures  1 through  3 show  the  steps 
taken  in  the  preparation  of  a silt  target. 

The  material  for  the  silt-cement  targets  was 
prepared  by  mixing  pre-wetted  silt  with  portland 
cement  in  a planetar)'  mixer.  The  amounts  of  port- 
land  cement  and  water  were  varied  to  yield  three 
different  formulations  having  about  the  same  con- 
sistency before  setting.  Water  and  cement  contents 
and  measured  strengths  for  the  three  formulations  arc 
given  in  Table  1. 

The  mixed  materials,  which  were  placed  in 
6-inch-diameter  by  12-inch-long  test  specimen  molds, 
were  cured  in  aspray-humidified  room  at  70°F.for  at 
least  2 weeks  prior  to  testing.  Strength  measurements 
were  obtained  within  a few  days  of  the  penetration 
tests.  The  cast-iron  molds  were  left  around  the 
specimens  during  the  penetration  testing  in  an 
attempt  to  prevent  the  specimens  from  splitting  so 
that  the  material  deformations  and  local  stresses  at 
the  penetrator  were  maintained  as  being  representa- 
tive of  the  behavior  in  a semi-infinite  space.  A silt- 
cement  target  after  testing  is  shown  in  Figure  4. 

Penetrators.  The  penetrators  used  in  the  silt 
model  tests  are  shown  in  Figures  5 and  6.  Half- 
section penetrators  (Figure  5)  were  1-3/8  inches 
in  diameter  by  approximately  12-1/2  inches  in 
length.  Two  half-section  steel  penetrators  were  used, 


one  with  a blunt  nose  and1  the  other  with  a conical 
nose  having  a length,  equal  to  twice  its  diameter.  A 
thin  layer  of  friction-reducing  plastic  material  was 
attached,  to  the  flat  side  of  each  half-section  pene- 
trator to  facilitate  movement  of  the  penetrator  along 
the  acrylic  plascic  face  of  the  test  tank  during  pene- 
tration. 

The  full-section  penetrators  used  in  the  silt  tests 
(Figure  6)  were  cylindrical  and  had  dimensions 
exactly  one-half  those  of  the  half-section  penetrators. 
A blunt-nosc  and  a pointed-nose  penetrator  were 
fabricated  from  both  steel  and  aluminum.  The  lighter 
material,  which  provides  greater  decelerations  under 
otherwise  similar  circumstances,  was  used  to  deter- 
mine the  effect  of  instantaneous  deceleration  on 
penetration  resistance. 

The  penetrators  used  in  the  silt-cement  and  con- 
crete model  penetration  tests  are  shown  in  Figures  7 
and  8.  The  cylindrical  penetrators  were  3/8-inch  in 
diameter  by  about  6 inches  in  length.  Conical  noses 
were  used  having  length-to-diameter  ratios  (sharp- 
nesses) of  2,  3,  and  4.  A “serrated”  quasi-cylindrical 
penetrator  also  was  used  having  a conical  nose  with  a 
sharpness  of  2,  followed  by  a body  consisting  of 
truncated  cones  each  having  a sharpness  of  4 and  a 
truncated  length  of  1/2  inch.  Also,  a lance-shaped 
penetrator  was  fabricated  having  a conical  nose  which 
tapered  from  a 1 /64-inch  diameter  at  the  point  to  a 
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1 /2-inch  diameter  at  'Ik-  rear  over  a length  of  about 
51/2  inches,  followed  by  a shaft  section  having  a 
3/8-inch  diameter.  The  Hat-plate  penctrators  (Figure 
8)  consisted  of  plates  l/I6-inch  thick  with  length-to- 
width  ratios  (slendernesses)  of  1,  1-1/2,  and  2;  each 
plate  had  a single-side  surface  area  of  approximately  6 
square  inches.  Two  additional  plates,  each  having  a 
thickness  of  1/16  inch,  a single-side  surface  area  of 
about  6 square  inches,  and  a slenderness  of  2,  were 
fashioned  into  a shape  similar  to  the  “sediment 
fluke"  anchor-projectile  of  an  existing  propcllant- 
h.unched  direct  embedment  anchor  (Taylor  anil 
Heard,  1973)  to  provide  data  on  the  performance  of 
that  anehor-piojectile.  All  of  the  projectiles  for  the 
cementc.l  materials  were  made  of  mild  steel,  ami  all 
had  a threaded  portion  approximately  1/2  inch  in 
length  at  the  rear  to  connect  it  to  a universal 
launching  shaft  and  instrument  housing  unit. 

Fquipmcnt.  The  model  test  tank  used  in  the  silt 
tests  (shown  in  Figure  3)  was  made  of  3/8-inch-thick, 
steel  painted  to  resist  corrosion.lt  consisted  of  a half- 
cylindrical  section  3 feet  in  diameter  by  about  3 feet 
in  length  attached  to  a half-hemispherical  end  section. 


Flanges  were  attached  to  all  edges  to  make  the  walls 
rigid  and  to  provide  a surface  for  bolting  and  sealing 
the  acrylic  plastic  face  to  the  half-sections.  A neo- 
prene rubber  seal  was  used  to  hold  the  water  that  was 
placed  several  inches  deep  over  the  top  of  each  silt 
target. 

The  penctrators  were  launched  from  a com- 
mercially available  stud  driver  gun  fitted  with  a barrel 
large  enough  to  drive  3/8-inch  studs.  The  gun  was 
modified  to  prevent  the  escape  of  the  propellant  gases 
lieforc  the  launching  shaft  had  cleared  the  barrel,  and 
a safety  feature  was  overridden  to  allow  the  gun  to  be 
fired  without  having  the  barrel  in  contact  with  any 
surface. 

The  penctrators  were  launched  downward 
through  a 3-inch-square  hole  cut  from  a 3/4-inch- 
thick  plywood  plate  placed  atop  the  test  tank.  The 
hole  allowed  some  lateral  adjustment  of  the  pene- 
trator  position  prior  to  firing,  and  it  provided 
sufficient  clearance  for  the  instrument  cable  to  follow 
the  penetrator  without  damage.  The  height  of  the  gun 
was  adjusted  with  plywood  spacers  so  that  the 
launching  shaft  would  exit  the  barrel  shortly  before 
the  penetrator  began  to  enter  the  target. 
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Figure  9.  Vacuum  pump  and  mixer  for  soil 
preparation. 


Several  pieces  of  equipment  were  required  for 
preparation  of  targets.  A planetary  mixer  was  used 
together  with  a bank  of  vacuum  pumps  as  shown  in 
Figure  9 to  provide  uniformly  mixed  and  de-aired  soil 
for  placement  in  the  tank.  Wooden  spacers  were  con- 
structed to  support  the  future  top  surface  of  the  tar- 
get as  soil  was  placed  in  the  tank  as  it  lay  on  its  side. 
A template  was  constructed  to  facilitate  placement  of 
talcum  powder  in  lines  1/8  inch  in  width  on  a 2 x 
2-inch  square  grid.  Other  equipment  utilized  included 
common  lifting  devices,  placement  tools,  containers, 
and  weighing  and  measuring  devices. 

Conventional  cast-iron  concrete  cylinder  molds 
were  utilized  for  the  cemented  targets.  The  cured 
targets  were  placed  underwater  in  the  silt  test  tank 
atop  a bed  of  compacted  sand.  After  testing,  an  air- 
driven  load-testing  machine  was  utilized  to. extract 
each  penetrator  from  its  target  while  the  extraction 
force  was  measured. 

Instrumentation.  High-speed  motion  pictures 
were  taken  with  a variable  speed  camara  at  rates  up  to 
5,000  frames  per  second.  Penetrator  travel  between 
frames  was  determined  by  comparing  the  movement 
of  the  projected  image  with  a scale  placed  on  the 


front  of  the  test  tank  which  also  appeared  on  the 
film;  overall  accuracy  was  approximately  0.1  inch  dif- 
ferential movement.  Timing  indicator  marks  were 
placed  on  the  film  at  0.001-second  intervals.  The 
operation  of  the  timing  indicator  marker  was  fairly 
good  during  the  early  test  scries,  but  it  deteriorated  in 
the  later  series;  it  produced  overall  penetrator  veloc- 
ity accuracies  varying  from  ±2%  to  as  much  as  ±30%. 

A piezoelectric  transducer  measured  the 
instantaneous  deceleration  of  the  penetrator  during 
each  test.  The  time  base  was  provided  by  the 
accurately  known  speeds  of  the  recording  devices. 
Three  different  models  of  the  accelerometer  were 
utilized  during  the  course  of  the  tests;  the  second  one 
was  substituted  for  the  first  one  to  reduce  size,  and 
the  third  one  replaced  the  second-ane  after  it  was 
damaged  by  excessive  shock  during  a test  into  silt- 
cement.  The  three  transducers  are  shown  in  Figure 
10.  The  third  one  (the  smallest)  proved  to  be  highly 
satisfactory  and  durable,  sustaining  peak  shock 
accelerations  in  excess  of  20,000  g.  One  of  the  trans- 
ducers is  shown  again  in  Figure  1 1 along  with  its 
mounting  screw  and  the  penetrator  adapter.  Also 
shown  are  the  launching  shaft,  accelerometer  housing, 
and  the  instrument  cable. 

The  output  of  the  accelerometer  was  fed 
through  a miniature  coaxial  cable  to  a charge 
amplifier.  The  cable  was  reasonably  durable  for  most 
tests,  but  behavior  became  erratic  for  the  tests  in  the 
concrete  targets.  The  charge  amplifier  conditioned 
the  signal  over  a wide  range  to  provide  a standardized, 
calibrated  output  to  a tape  recorder.  A parallel  signal 
from  the  charge  amplifier  was  filtered  and  monitored 
on  an  oscilloscope,  from  which  a permanent  record 
was  obtained  with  a still-  camera  mounted  on  the 
oscilloscope  screen. 

The  rctual  distances  traveled  by  the  models  were 
measured  with  a steel  tape  accurate  to  within  approx- 
imately 0.03  inch  over  the  distances  measured. 

A final  instrumentation  item  was  the 
load-testing  machine  utilized  to  extract  the  pene- 
trators  from  the  cemented  targets.  Kxtraction  forces 
usually  were  monitored  on  the  low  range  of  this 
machine  (up  to  600  pounds  full  scale  with  an 
accuracy  of  approximately  ±10  pounds).  When  the 
extraction  load  exceeded  600  pounds,  the 
6,000-pound,  full-scale  range  was  utilized  with  an 
accuracy  on  the  order  of  ±50  pounds. 
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Data  Reduction 

Two  types  of  measurement  required  reduction 
to  a more  meaningful  form:  the  high-speed  motion 
pictures  and  the  accelerometer  data.  The  high-speed 
movies,  in  addition  to  providing  qualitative  informa- 
tion on  soil  deformation  patterns,  enabled  penettator 
velocity  to  be  determined  by  scaling  distances  and 
calculating  film  speed.  Film  speed  was  determined 
from  the  measured  distances  between  the  timing 
indicator  marks.  The  timing  marks  (recorded  on  the 
film  at  I /msec)  were  numbered,  beginning  with  zero 
at  the  stan  of  noticeable  pcnctrator  motion.  The 
distances  to  the  timing  marks  were  scaled  in  terms  of 
number  (and  decimal  fraction)  of  frames.  The  slope 
for  a plot  of  number  of  frames  versus  timing  mark 
number  was  used  to  compute  the  film  speed  in  frames 
per  second.  Then,  from  the  computed  film  speed, 
precise  times  were  computed  for  each  film  frame.  The 
scaled  position  of  the  penetrator  was  plotted  against 
this  computed  time  for  each  film  frame.  The  pene- 
trator velocity  was  determined  from  the  slope  of  this 
plot,  with  the  impact  velocity  corresponding  to  the 
greatest  slope. 

The  accelerometer  data,  recorded  as  an  analog 
frequency  on  magnetic  tape,  was  digitized,  input  into 
a computer,  and  integrated  twice  with  respect  to  time 
to  obtain  velocity  and  displacement.  Velocity  was 
augmented  by  a constant  so  that  it  diminished  to  zero 
at  the  end  of  the  event,  signified  by  the  acceleration 
trace  ceasing  to  show  deviation  from  a constant 
(zero)  value.  Displacement  was  taken  as  zero  at  the 
start  of  penetration,  which  was  readily  detectable 
from  the  acceleration  trace.  Deceleration,  velocity, 
and  displacement  were  plotted  against  time  and 
against  each  other,  two  at  a time,  to  give  a total  of  six 
plots.  These  plots  were  used  to  determine  values  of 
instantaneous  (as  opposed  to  average)  deceleration 
versus  velocity  and  displacement  for  several  equally 
spaced  points  during  each  event;  this  information  was 
further  analyzed  by  hand. 

Although  precautions  were  taken  to  ensure  that 
all  instrumentation  would  function  properly,  the 
performance  of  the  acceleration-measuring  instrumen- 
tation on  some  tests  required  the  data  to  be  modified 
to  fit  other  measurements.  The  problem  involved 
drift  of  the  zero-acceleration  output  of  the  system 
with  time,  especially  when  the  charge  amplifier  was 


set  on  a sensitive  range  (as  for  the  s3t  tests).  This 
problem,  which  often  occurs  with  piezoelectric  trans- 
ducers, was  caused  by  one  or  more  small  electrical 
leaks  ip  the  circuits  upstream  from  the  output  of  the 
chaige  amplifier.  If  leaks  occur  in  the  acceleration- 
sen skive  dements,  such  as  the  cable  or  transducer,  a 
shift  in  zero  output  can  occur  abruptly  when  large 
accelerations  are  sustained  by  these  elements. On  the 
other  hand,  leaks  in  elements  not  sensitive  to  acceler- 
ation produce  a slow  drift.  Smooth  drifting  was 
observed  during  some  of  the  silt  tests,  and  abrupt 
shifts  were  observed  during  some  of  the  cemcnted- 
material  tests.  The  cables  were  replaced  and 
precautions  were  taken  to  reduce  the  likelihood  of 
obtaining  erratic  data,  with  some  success.  However, 
successful  data  interpretation  required  correction  of 
the  erratic  data  by  adjustment  of  zero  levels  in  the 
computer  operations.  Further,  when  the  computed 
overall  velocity  and  displacement  changes  were 
checked  against  film-measured  velocities  and 
manually  measured  displacements,  the  latter  were 
given  more  weight  and  were  used  to  correct  computer 
values  by  a proportionate  adjustment  of  instan- 
taneous data. 

Test  Data 

A typical  view  of  a silt  target  after  penetration 
by  a half-section  projectile  is  shown  in  Figure  12.  The 
soil  deformation  patterns  that  were  observed  through- 
out motion  pictures  of  such  penetration  events  arc 
shown  in  Figure  13.  Also,  high-speed  films  provided 
penetrator-position-versus-time  data;  these  data  were 
plotted,  and  the  slopes  of  the  data  curves  were  used 
to  compute  impact  velocity  for  correcting  the 
accelerometer  data. 

A deceleration  trace  from  the  oscilloscope 
screen  is  shown  in  Figure  14  for  a silt  test  with  a 
blunt-nosed  aluminum  pcnctrator.  The  computer- 
reduced  plots  for  this  test  arc  shown  in  Figure  15. 

The  test  data  from  all  of  the  model  tests  arc 
summarized  in  Tabic  2 along  with  analytical  calcula- 
tions that  will  be  discussed  in  a following  section. 
These  data  arc  overall-type  parameters,  such  as 
impact  velocity  and  ultimate  depth  of  embedment. 

Instantaneous  conditions  were  evaluated  for 
several  points  in  the  soft  and  medium  silt  and  in  the 
silt-cement  test  scries.  These  data  arc  summarized  in 
Tables  3,  4,  and  5,  along  with  analytical  computa- 
tions that  will  be  discussed  in  a following  section. 
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Figure  13.  Artist’s  reconstruction  of  soil  deformations  observed  in  high-speed 
movies;  silt  target  with  blunt  half-section  pcnctrator. 


Figure  14.  Deceleration  trace  of  a blunt-nosed  aluminum  pcnctrator  in 
silt-target  test. 
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Symbol 

Fenctrator 

Used  in 

Designation 

Plotring 

Data 

pointed 

• 

pointed 

• 

pointed 

• 

pointed 

pointed 

blunt 

blunt 

blunt 

blunt 

blunt  steel 
pointed  aluminum 
blunt  aluminum 
blunt  steel 
blunt  steel 
pointed  steel 
pointed  aluminum 
pointed  steel 

blunt  steel 
pointed  steel 
blunt  aluminum 
pointed  aluminum 
blunt  steel 
pointed  steel 


pointed 

pointed  serrated 
pointed 
very  pointed 
lance 

2:1  plate 
MAM  plate 
1:1  plate 
pointed 
pointed 
very  pointed 
lance 


Effective  Areas 


Side  wall  serrated. 


m 

1SK71 


Material 

Sensitivity, 

St 


Material 
Density, 
7 » Pf 
(gf/cm*) 


Bearing 

Pressure 

Factor 


LSP1 

very  pointed 

I.SP3 

pointed 

1,250,000 

0.31 

1,250,000 

0.14 

2.39 

80.3 

2.39 

80.3 

2.39 

80.3 

2.39 

80.3 

2.39 

80.3 

2.39 

75.4 

2.39 

75.4 

2.39 

75.4 

2.39 

75.4 

2.39 

68 

2.18 

59 

2.39 

59 

2.39 

65 

2.39 

61 

2.15 

55 

2.16 

59 

2.23 

63 

2.39 

67 

2.18 

54 

2.39 

58 

2.17 

51 

2.39 

64 

2.16 

48 

1.4 

75 

0.9 

46 

0.8 

30 

>1.6 

>45 

0.70 

23.5 

0.65 

11.5 

0.6 

18.5 

0.57 

14 

0.7 

10.3 

0.8 

8.3 

5.0 

8.9 

0.40 

2.5 

0.52 

5.0 

0.20 

1.0 

0.18 

4.0 

1.6 

11 

170 

1.5 

11 

170 

v = v0  for  static  tests,  2/3  v0  for  dynamic  tests. 


Table  2.  Overall  Test  Data  and  Computations. 


ectiv  Areas 
(an2) 

Effective 

Maas, 

an 

Effective 

Velocity. 

(m/sec) 

Embedment 
Depth, Z 

(cm) 

— 

Effective 
Length  of 
Shear  Zooe, 
L(cm) 

Effective 

Avenge 

Decektatioo. 

“a 

(g)Z 

Penetration 

Rcaauace 

Force, 

F > M*i  + W ♦ Fp 
Otgf) 

Net  Soil 
Resaance, 
F*  F FH 

dp  dp  dp 

(tgf/cm ?) 

Normalized 

Velocity, 

* 

i 

Side, 

(10'2  m2/kgf-aec) 

Laboratory  Model  Teats 

9 

80.3 

■ 

0.232 

Static 

3.00 

1.259 

0.862 

9 

80.3 

- 

0.219 

- 

Sutic 

2.64 

1.109 

0.811 

9 

80.3 

— 

1.140 

— 

15.6 

Sutic 

3.36 

1.391 

4322 

9 

80.3 

- 

0.960 

— 

15.0 

Static 

3.77 

1.567 

3JS 

9 

80.3 

— 

031257 

— 

15.0 

Static 

1.68 

0.705 

03)944 

9 

75.4 

- 

0.0265 

- 

15.0 

Static 

1.24 

0.519 

03)987 

9 

75.4 

- 

0.0656 

- 

153) 

Sutic 

2.88 

1.199 

2.43 

9 

75.4 

- 

0.254 

— 

153) 

Static 

3318 

1.289 

0.938 

9 

75.4 

- 

0.272 

- 

15.0 

Sutic 

2.16 

0.904 

1.004 

9 

68 

431 

633 

68.9 

13.4 

3.47 

1.84 

0.635 

18.10 

8 

59 

263 

8.75 

37.5 

12.0 

11.15 

3.11 

1.186 

26.75 

9 

59 

262 

7.20 

45.5 

12.5 

5.78 

1.68 

0.541 

21.12 

9 

65 

431 

5.35 

49.0 

12.7 

3.63 

1.91 

0.707 

15.51 

9 

61 

262 

6j62 

37.8 

12.0 

6.42 

1.85 

0.635 

20.36 

5 

55 

434 

6.54 

29.9 

11.2 

7.40 

3.55 

1.514 

21.39 

6 

59 

263 

9.08 

33.7 

11.7 

12.50 

3.46 

1.341 

28.71 

3 

63 

434 

8.23 

41.3 

12.3 

8.65 

4.10 

1.624 

24.65 

9 

67 

438 

>12.43 

>73.8 

>13.5 

<11.5 

<5.37 

<1.725 

>17.69 

8 

54 

441 

8.14 

30.0 

11.3 

12.6 

5.90 

2.485 

13.76 

9 

58 

269 

9.20 

33.9 

11.7 

14.2 

3.99 

1.382 

14.98 

7 

51 

270 

12.95 

26.6 

10.8 

34.7 

9.54 

33330 

17.42 

9 

64 

438 

8.22 

38.1 

12.0 

10.6 

4.98 

1.850 

13.02 

6 

48 

441 

9.04 

22.9 

10.1 

19.2 

8.82 

3.802 

16.88 

75 

309 

49.8 

23.13 

6.62 

545 

169 

111 

0.0421 

46 

233 

50.8 

10.93 

5.07 

1,201 

280 

302 

0.0124 

30 

222 

51.8 

7.11 

3.50 

1.922 

425 

521 

0.00449 

>45 

>351 

51.5 

>30.7 

<11.34 

440 

157 

<88 

0.0253 

m 

23.5 

222 

53.4 

17.78 

8.66 

817 

182 

249 

0.00762 

5 

11.5 

214 

50.8 

9.65 

3.34 

1,361 

291 

437 

0.00461 

18.5 

219 

53.1 

16.25 

8.05 

884 

194 

311 

0.00200 

14 

221 

54.4 

15.5 

7.75 

973 

215 

365 

0.00212 

10.3 

201 

54.8 

2.24 

1.12 

6,830 

1,372 

1,950 

0.00808 

8.3 

197 

50.3 

1.58 

0.79 

8,150 

1,607 

1,995 

0.00480 

8.9 

204 

57.4 

1.37 

0.69 

12,230 

2,500 

2,488 

0.00324 

2.5 

196 

23.3 

2.92 

1.46 

948 

186 

463 

0.00263 

5.0 

198 

67.2 

4.7 

2.35 

4,890 

968 

1.845 

0.00472 

1.0 

187 

61.7 

3.18 

1.59 

6,100 

1,140 

5,686 

0.00150 

4.0 

196 

69.6 

7.15 

3.58 

3,450 

676 

3,730 

0.00146 

Large  Seafloor  Penetrometer  Tests 

m 

1,562,000 

17 

■m 

400 

3.95 

7,260 

■O 

1 

H 

1,543,000 

7 

350 

0.88 

2,454 

> 

i 

•s 

Rniifft 

Force, 

» uH*w*rn 


Net  Soil 
Resistance, 

F«  F FH 

% *r  *f 

(k*f/cm2) 


Minimum  Strain -Rate 
Effect. 

FVAp 


Vc^VtXVV 


Beck -Calculated 
Side  Adhesion 
Factor, 6 


3.00 

1.259 

0362 

6.35 

2.09 

0.790 

2.64 

1.109 

0311 

5.60 

134 

0.654 

3.36 

1.391 

432 

7.03 

2.31 

0.534 

3.77 

1.567 

335 

7.92 

2.60 

0.69$ 

1.6S 

0.705 

0.0944 

3.56 

1.17 

0.809 

1.24 

0319 

0.0987 

3.20 

0.96 

0.582 

22t 

1.199 

2.43 

7.40 

2.21 

0.606 

3j08 

1.289 

0.938 

7.95 

2.38 

0.940 

2.16 

0.904 

1.004 

538 

1.67 

0.516 

1.84 

0.635 

18.10 

3.92 

1.35 

-0.017 

3.11 

1.186 

26.75 

5.99 

2.14 

0.239 

1.6B 

0341 

21.12 

3.34 

1.26 

-0.112 

1.91 

0.707 

15.51 

4.37 

1.55 

0.052 

1.85 

0.635 

20.36 

3.92 

1.45 

-0.025 

335 

1314 

21.39 

7.64 

2.83 

0.505 

3.46 

1.341 

28.71 

6.76 

2.40 

0.333 

4.10 

1.624 

24.65 

8.20 

2.85 

0.516 

<5.37 

<1.725 

>17.69 

5.47 

<1.91 

<0.189 

5.90 

2.485 

13.76 

6.45 

2.39 

0.376 

3.99 

1.382 

14.98 

4.J9 

1.68 

0.064 

9.54 

3.830 

17.42 

10.00 

3.74 

0.885 

4.98 

1.850 

13.02 

5.89 

2.12 

0.282 

8.82 

3.802 

16.88 

9.90 

3.98 

0.984 

169 

111 

0.0421 

13.3 

2.92 

0.829 

280 

302 

0.0124 

8.0 

2.95 

1.178 

425 

521 

0.00449 

3.36 

1.96 

0.886 

157 

<88 

0.0253 

6.33 

<2.62 

0.748 

182 

249 

0.00762 

4.19 

2.69 

1.386 

291 

437 

0.00461 

1.79 

1.43 

194 

311 

0.00200 

1.17 

0.62 

215 

365 

0.00212 

1.10 

0.36 

1,372 

1,950 

0.00808 

6.90 

4.05 

2.668 

1.607 

1.995 

0.00480 

3.25 

2.17 

1.115 

2,500 

2.488 

0.00324 

2.07 

1.45 

0.236 

186 

463 

0.00263 

1.04 

0.87 

-2.177 

968 

1.845 

0.00472 

4.15 

3.22 

3.407 

1,140 

5,686 

0.00150 

2.84 

2.49 

6.141 

676 

3.730 

0.00146 

2.83 

1.90 

1.775 

Deceleration  (g)  Deceleration  (g) 
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Figure  1 5.  Computer-processed  data  for  a typical  silt-target  model  test. 
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sc  Data  and  Computations  for  Models  in  Soft  Silt. 


. Normalized 
Velocity, 

V.*JL 

LS- 

(m2/kgf-sec) 

Fluid 

Inertial  Drag, 

FH  *2^2cDAF 
(kgf) 

Net  Measured 
Penetration 
Resistance, 
F*=M*a-FH  +W* 
(kgf) 

Maximum 

Strain-Rate 

Effect, 

fvaf 

S*«ax  = S„NC 

Minimum 
Strain-Rate  Effect, 
F*/Ap 

«mm  SBNc  + (Su/St)(As/AF) 

0.246 

0.45 

1.38 

3.56 

1.02 

0.188 

0.31 

1.59 

4.11 

1.17 

0.151 

0.20 

1.40 

3.62 

1.03 

6.113 

0.11 

1.59 

4.11 

1.17 

0.075 

0.05 

2.02 

5.22 

1.49 

0:038 

0.01 

2.15 

5.55 

1.59 

0.429 

0.83 

1.09 

2.30 

0.97 

0.277 

0.64 

3.36 

7.10 

2.23 

0.238 

0.47 

3.47 

7.33 

2.40 

0.198 

0.32 

3.56 

7.51 

2.36 

0.158 

0.21 

3.62 

7.65 

2.41 

0.119 

0.12 

3.60 

7.61 

2.39 

0.079 

0.05 

3.56 

7.52 

2.37 

0.040 

0.01 

3.23 

6.83 

2.14 

0.240 

0.49 

1.48 

3.82 

1.09 

0.199 

0.34 

1.43 

3.70 

1.06 

0.159 

0.22 

1.61 

4.16 

1.19 

0.119 

0.12 

1.68 

4.34 

1.24 

0.079 

0.06 

1.72 

4.44 

1.27 

0.040 

0.01 

1.93 

4.99 

1.42 

1.643 

0.34 

0.80 

2.06 

1.54 

0.158 

0.22 

1.06 

2.74 

0.78 

0.118 

0.12 

2.40 

6.20 

1.77 

0.079 

0.05 

3.07 

7.94 

2.27 

0.039 

0.01 

3.02 

7.80 

2.23 

0.602 

0.51 

-0.02 

•0.05 

-0.02 

0.201 

0.35 

1.30 

3.36 

0.96 

0.161 

0.23 

1.65 

4.26 

1.22 

0.121 

0.13 

2.27 

5.86 

1.68 

0.080 

0.06 

2.46 

6.35 

1.82 

0.040 

0.01 

2.53 

6.54 

1.87 

continued 


Table  3.  Continued 


3.  Continued 


lormalizcd 

Velocity, 

LSu 

i2/kgf-sec) 

Fluid 

Inertial  Drag, 

fH  = 2°v2cdaF 
(kgO 

Net  Measured 
Penetration 
Resistance, 

f*=m*«-fh  + w' 

(kgf) 

Maximum 

Strain-Rate 

Effect, 

F*/Af 
s*  = F 

«mix  SUNC 

. 

Minimum 
Strain-Rate  Effect, 

f*/af 

“emin  SUNC  + <Su/Sc)(As/Af) 

0.343 

0.44 

2.52 

5.33 

2.37 

0.179 

0.28 

3.74 

7.91 

2.49 

0.135 

0.16 

4.27 

9.04 

2.84 

6.090 

0.07 

4.55 

9.62 

3.03 

0.045 

0.02 

4.12 

8.72 

2.74 

0.440 

0.87 

1.75 

3.70 

1.42 

0.277 

0.67 

3.22 

6.80 

2.13 

0.236 

0.49 

3.97 

8.40 

2.64 

0.198 

0.34 

4.55 

9.62 

3.02 

0.158 

0.22 

4.67 

9.88 

3.11 

0.119 

0.12 

4.31 

9.12 

2.87 

0:079 

0.05 

4.10 

8.67 

2.73 

0.040 

0.01 

3.44 

7.27 

2.28 

0.263 

0.60 

3.54 

7.49 

2.37 

0.210 

0.38 

4.39 

9.28 

2.98 

0.158 

0.22 

4.83 

10.21 

3.21 

0.105 

0.10 

5.08 

10.74 

3.38 

0.053 

0.02 

4.56 

9.65 

3.04 

a and  Computations  for  Models  in  Medium  Silt. 


formalized 

Velocity, 

"** 

n2/kgf-sec) 

Fluid 

Inertial  Drag, 
(kgO 

Net  Measured 
Penetration 
Resistance, 

F*  sM’s-Fm  + W1 

<k*0 

Maximum 

Strain'Rate 

Effect, 

F*/AF 

S‘u«  ' S„NC 

Minimum 
Strain*Rate  Effect, 

F*/Af 

s-  - 

S,Nc+(S./StKAs/AF) 

0.765 

1.90 

2.06 

2.74 

1.57 

0.228 

1.60 

3.27 

4.35 

1.60 

0.208 

1.32 

4.82 

6.40 

1.84 

0.187 

1.07 

5.35 

7.11 

2.10 

0.166 

0.84 

5.76 

7.65 

2.26 

0.146 

0.65 

6.08 

8.08 

2.38 

0.125 

0.47 

6.03 

8.01 

2.36 

0.104 

0.33 

5.90 

7.84 

2.31 

0.083 

0.21 

6.48 

8.60 

2.54 

0.062 

0.12 

7.29 

9.68 

2.86 

0.042 

0.05 

6.64 

8.82 

2.60 

0.021 

0.01 

5.13 

6.81 

2.01 

0.191 

0.68 

4.12 

4.48 

1.55 

0.127 

0.55 

6.05 

6.58 

2.24 

0.113 

0.44 

6.92 

7.53 

2.48 

0.099 

0.33 

7.30 

7.94 

2.61 

0.085 

0.25 

7.51 

8.16 

2.69 

0.070 

0.17 

7.76 

8.44 

2.78 

0.056 

0.11 

7.85 

8.53 

2.81 

0.042 

0.06 

7.63 

8.30 

2.97 

0.028 

0.03 

6.96 

7.57 

2.71 

0:014 

0.01 

6.25 

6.80 

2.43 

1.52 

1.05 

1.90 

2.52 

1.95 

0.271 

0.85 

2.87 

3.81 

1.55 

0.145 

0.67 

2.68 

3.57 

1.03 

0.121 

0.51 

2.66 

3.53 

1.03 

0.104 

0.38 

3.34 

4.44 

1.29 

0.086 

0.26 

5.21 

6.93 

2.01 

0.069 

0.17 

5.79 

7.69 

2.27 

0.052 

0.09 

5.99 

7.95 

2.31 

0.035 

0.04 

6.09 

8.09 

2.35 

0.017 

0.01 

5.31 

7.05 

2.05 

0.404 

1.88 

4.65 

5.05 

2.55 

0.270 

1.62 

7.89 

8.68 

3.53 

0.198 

1.38 

10.19 

11.07 

3.97 

0.182 

1.16 

11.11 

12.10 

3.98 

0.165 

0.96 

11.56 

12.57 

4.14 

continued 
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Table  4.  Continued 


Tabic  4.  Continued 


Side  Area, 

AS 

(cm2) 


Normalized 
Velocity, 
v 


v*  = 


80.3 

80.3 

80.3 

80.3 

80.3 

80.3 

80.3 

80.3 

80.3 

4.5 

71.7 

75.4 
75.4 
75.4 
75.4 
75.4 
75.4 
75.4 

75.4 

2 

35 

53 

67.5 
79 
80.3 
80.3 
80.3 
80.3 
80.3 
80.3 


LS. 


(m2/kgf-*c) 


0.149 

0.132 

0.115 

0.099 

0.082 

0.066 

0.049 

0.033 

0.017 

3-03 

0.160 

0.125 

0.110 

0.094 

0.078 

0.063 

0.047 

0.031 

0.016 

1.709 

0.295 

0.185 

0.134 

0.109 

0.093 

0.078 

0.062 

0.047 

0.031 

0.016 


Fluid 

Inertial  Drag, 

fh  =7^’cdaf 

0*0 


0.78 

0.61 

0.47 

0.35 

0.24 

0.15 

0.09 

0.04 

0.01 

0.85 

0.69 

0.55 

0.42 

0.31 

0.21 

0.14 

0.08 

0.03 

0.01 

1.03 

0.85 

0.69 

0.54 

0.42 

0.31 

0.21 

0.14 

0.08 

0.03 

0.01 


Net  Measured 
Penetration 
Resistance, 
F*=M*a-F„+W' 
(kgf) 


12.02 

12.43 

12.63 

12.64 
12.33 
12.09 
11.71 
10.98 

9.16 

0.91 

1.73 

2.41 

5.52 

6.51 

7.19 
7.48 
7.83 
7.83 
7.43 

1.88 

6.20 
8.12 
9.19 
9.92 

10.37 

10.43 
10.43 

10.16 
9.70 
8.58 


Maximum 

Strain-Rate 

Effect. 

¥m/Ar 


V»c 


13.08 

13.52 

13.74 

13.75 

13.40 
13.14 
12.74 
11.93 

9.96 

1.21 

2.30 

3.20 

7.34 

8.65 

9.55 

9.93 

10.40 
10.47 

9.86 

12.20 

6.74 

8.83 

9.99 

10.80 

11.27 

11.35 

11.35 

11.03 

10.54 

9.32 


Minimum 
Strain-Rate  Effect, 

FVAf 

= S.Nc+(S>/StKAs/AF) 


4.30 

4.45 

4.52 

4.53 
4.42 
4.33 
4.20 

3.93 
3.28 

1.06 

0.72 

0.93 

2.13 

2.52 

2.78 

2.89 

3.02 

3.04 

2.87 

1.94 
3.57 
3.76 
3.68 
3.59 
3.72 
3.74 
3.74 
3.64 
3.48 
3.07 


Tabic  5.  Instantaneous  Test  Data  and  Computations  for  Models  in  Silt 


Velocity, 

v 

(cm/sec) 


Deceleration, 

a 

<g> 


Penetration 
Depth,  z 
(cm) 


Frontal  Area, 

af 

(cm2) 


Side  Area, 
Ac 


Fluid 

Inertial  Drag, 

fh  *‘JPv2cdaf 
(kgO 


8.48 

0.925 

75 

— 

23.3 

13.00 

1.20 

86.5 

....  1 

22.2 

16.67 

1.20 

93 

0.0277 

15.45 

19.4 

3.255 

‘ 107.5 

>0.0186 

26.8 

20.9 

3.255 

115.8 

>0.0129 

15.0 

22.1 

3.255 

123 

>0.0082 

6.65 

22.8 

3.255 

125 

>0.0040 

1.66 

2.96 

1 

0.705 

23 

0.0398 

16.2 

4.84 

0.925 

41 

0.0224 

17.2 

6.92 

0.925 

60 

0.0131 

12.5 

8.44 

0.925 

75 

0.0086 

8.0 

9.49 

1.20 

75 

0.0057 

5.9 

10.22 

1.20 

81 

0.0036 

2.6 

10.50 

1.20 

81.4 

0.0026 

1.5 

10.65 

1.20 

81.7 

0.0011 

0.7 

10.80 

1.20 

82 

0.9008 

0.2 

i 

2.07 

0.705 

15 

0.0109 

10.4 

3.32 

0.925 

26 

0.0061 

11.1 

4.17 

0.925 

29.5 

0.0043 

8.7 

4.95 

0.925 

42 

0.0031 

6.7 

5.55 

0.925 

48 

0.0024 

4.9 

6.02 

0.925 

52 

0.0019 

3.4 

6.37 

0.925 

55 

0.0014 

2.2 

6.60 

0.925 

57.5 

0.0010 

1.2 

6.80 

0.925 

59 

0.0007 

0.5 

6.93 

0.925 

60.5 

0.0003 

0.1 

3.24 

0.713 

7 

0.228 

32.7 

15.84 

0.713 

44 

0.0415 

25.8 

22.1 

2.852 

76.7 

0.0260 

79.1 

25.1 

2.852 

76.7 

0.0196 

58.1 

27.7 

2.852 

76.7 

0.0148 

40.4 

30.2 

2.852 

76.7 

0,0109 

25.8 

Data  and  Computations  for  Modcb  in  Silt -Cement 


Normalized 

Velocity, 

v*«-L 

(m2/kgf-stc) 

Fluid 

Inertial  D>*g, 

fh  x7'*r2cDAF 
(kgO 

Net  Measured 
Penetration 
Resistance, 
F*«M*a-FM-iV 
0*0 

Maximum 

Strain-Rate 

Effect, 

f*/af 

$,NC 

Minimum 
Strain-Rate  Effect. 
FVAf 

«.  r 

c»«  S.Nc  + (S||,'S|KAs/Af) 

0.0745 

23.3 

181 

21.6 

3.18 

0.0416 

212 

287 

28.6 

4.66 

0.0277 

15.45 

343 

34.2 

5.30 

>0.0186 

26.8 

378 

14.0 

3.91 

>0.0129 

15.0 

418 

15.5 

4.38 

>0.0082 

6.65 

392 

14.5 

3.92 

>0.0040 

1.66 

366 

13.6 

3.64 

0.0398 

16.2 

182 

6.81 

3.14 

0.0224 

17.2 

297 

8.47 

3.28 

0.0131 

12.5 

396 

11.3 

3.40 

1 0.0086 

8.0 

458 

13.1 

3.37 

! 0.0057 

5.9 

472 

10.4 

3.10 

0.0036 

2.6 

462 

10.2 

2.99 

0.0026 

1.5 

445 

9.78 

2.86 

0.0011 

0.7 

388 

8.53 

2.48 

0.0008 

0.2 

286 

6.29 

1.83 

0.0109 

10.4 

287 

2.71 

1.89 

0.0061 

l'U 

419 

2.91 

1.85 

0.0043 

8.7 

424 

2.95 

1.79 

0.0031 

6.7 

502 

3.49 

1.81 

0.0024 

4.9 

559 

3.88 

1.88 

0.0019 

3.4 

612 

4.26 

1.98 

0.0014 

2.2 

651 

4.53 

2.04 

0.0010 

1.2 

645 

4.49 

1.98 

0.0007 

0.5 

620 

4.31 

1.87 

0.0003 

0.1 

589 

4.10 

1.76 

0.228 

32.7 

57 

6.09 

4.21 

0.0415 

25.8 

224 

23.9 

6.28 

0.0260 

79.1 

349 

9.32 

4.19 

0.0196 

58.1 

513 

13.7 

6.16 

0.0148 

40.4 

495 

13.22 

5.95 

0.0109 

25.8 

206 

5.50 

2.47 

continued 
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Tabic  5.  Continued 


Test  No. 

Symbol 
Used  in 
Plotting 
Dan 

Velocity, 

V 

(cm/sec) 

Deceleration, 

a 

<g> 

Penetration 
Depth,  a 
(cm) 

. 

Frontal  Area, 

V 

(an2) 

Side  Area, 
(an2) 

Normalized 

Velocity, 

(m2/kgf-scc) 

Fkid 

Inertial  Drag, 
FH 

(kgf) 

R3B 

4.900 

4.63 

— 

10 

0.0215 

11.7 

4,400 

7.8 

19.5 

9.7 

3,570 

1,320 

11.5 

0.713 

26 

0.0063 

62 

2,670 

1,350 

14.1 

0.713 

38.5 

35 

1,780 

1,380 

16.1 

0.713 

44.5 

1 

1.5 

890 

1,160 

17.2 

0.713 

47.5 

0.4 

450 

840 

17.6 

0.713 

48.5 

0.0005 

0.1 

R4B 

• 

4,770 

1,120 

2.92 

0.713 

5 

0.0076 

10.5 

4,290 

1,580 

4.42 

0.713 

9 

0.0045 

8.5 

3,820 

1,720 

5.60 

0.713 

13 

6.7 

3,340 

1,780 

6.65 

0.713 

16.5 

0.0023 

52 

2,860 

1,800 

7.36 

0.713 

18.5 

32 

2,380 

1,670 

8.06 

0.713 

20.5 

0.0014 

2.6 

1,910 

1,640 

8.62 

0.713 

22.5 

0.0010 

1.7 

1,430 

1,760 

9.10 

0.713 

0.0007 

1.0 

950 

1,750 

9.36 

0.413 

24 

0.0005 

0.4 

480 

1,590 

9.56 

0.713 

24.5 

0.0002 

0.1 

R5B 

• 

5,050 

1,150 

1.91 

0.75 

0.0116 

3-3 

4,540 

960 

4.76 

7 

0.0042 

8.5 

4,040 

920 

7.76 

0.713 

16 

0.0023 

6.7  • 

3,530 

1,080 

10.06 

0.713 

22.5 

0.0015 

5.1 

3,030 

1,130 

11.54 

0.713 

27 

0.0012 

3.8 

2,520 

1,250 

13.3 

0.713 

32 

0.0008 

2.6 

2,020 

1,260 

14.3 

0.713 

35 

0.0006 

1.7 

1,510 

1,230 

15.2 

0.713 

37.7 

0.0004 

0.9 

1,010 

1,530 

15.9 

0.713 

40 

0.0003 

0.4 

500 

1,610 

16.2 

0.713 

40.7 

0.0001 

0.1 

R5A 

A 

5,350 

350 

1.61 

0.061 

0.57 

0.0124 

0.7 

4,760 

800 

6.45 

0.325 

7.61 

0.0027 

3.1 

4,170 

1,030 

9.20 

0.570 

14.1 

0.0017 

4.1 

3,570 

1,200 

11.20 

0.790 

19.9 

0.0012 

4.1 

2,980 

1,350 

12.57 

0.960 

24.5 

0.0009 

3.5 

2,380 

1,430 

13.7 

1.115 

28.6 

0.0006 

2.6 

1,780 

1,500 

14.46 

1.205 

31 

0.0005 

1.6 

1,190 

1,700 

15.1 

1.205 

31 

0.0003 

0.7 

590 

1,690 

15.36 

1.205 

31 

0.0001 

0.4 

TableS.  Continued 


Side  Area, 
(«£) 

Normalized 

Velocity, 

v.-JL 

LS. 

(m2/kgf-*c) 

Fluid 

Inertial  Draj, 
fh  “ 

(kgf) 

Net  Measured 
Penetration 
Resistance, 

F*  = M!a-F„+W* 
<1*0 

Maximum 

Strain-Rate 

Effect, 

F*/Af 

S>NC 

Minimum 
Strain-Rate  Effect, 
FVA, 

s*«-  * s,n£ + <s„/stx VV 

10 

0.0215 

11.7 

211 

■I 

?.76 

19.5 

9.7 

257 

3.73 

26 

0.0063 

62 

287 

3.69 

38.5 

0.0039 

3.5 

297 

7.00 

3.14 

44.5 

1.5 

305 

7.18 

2.97 

47.5 

0.4 

257 

6.05 

2.41 

48.5 

0.0005 

0.1 

187 

4.40 

1.73 

5 

0.0076 

10.5 

229 

1.32 

1.21 

9 

0.0045 

8.5 

330 

1.90 

1.63 

13 

0.0032 

6.7 

362 

2.08 

1.68 

16.5 

0.0023 

5.2 

376 

2.16 

1.66 

18.5 

0.0018 

3.8 

381 

2.19 

1.64 

20.5 

0.0014 

2.6 

355 

2.04 

1.49 

22.5 

0.0010 

1.7 

350 

2.01 

1.61 

0.0007 

1.0 

375 

2.15 

1.51 

24 

0.0005 

0.4 

375 

2.15 

1.50 

24.5 

0.0002 

0.1 

340 

1.95 

1.35 

0.75 

0.0116 

3.3 

249 

4.17 

4.01 

7 

0.0042 

8.5 

202 

1.06 

0.95 

16 

6.7 

195 

1.03 

0.81 

22.5 

5.1 

232 

1.22 

0.89 

27 

0.0012 

3.8 

244 

1.28 

0.89 

32 

0.0008 

2.6 

272 

1.43 

0.93 

35 

0.0006 

1.7 

275 

1.45 

0.91 

37.7 

0.0004 

268 

1.41 

0.87 

40 

0.0003 

0.4 

335 

1.76 

1.06 

40.7 

0.0001 

0.1 

353 

1.85 

1.11 

0.57 

0.7 

77 

3.80 

3.49 

7.61 

3.1 

174 

1.62 

1.32 

14.1 

4.1 

224 

1.24 

1.00 

19.9 

0.0012 

4.1 

261 

1.00 

0.81 

24.5 

0.0009 

3.5 

295 

0.93 

0.75 

28.6 

0.0006 

2.6 

314 

0.85 

0.68 

31 

0.0005 

1.6 

330 

0.83 

0.66 

31 

0.0003 

0.7 

375 

0.94 

0.75 

31 

0.0001 

0.4 

373 

0.93 

0.74 
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FIELD  TESTS 


Fice-Fal  Coae  wd  Cwer 

Tests  were  conducted  early  in  the  course  of  the 
CEL  penetration  work  with  a free-fall  device  instru- 
mented to  record  deceleration  versus  time  on  a 
rotating  drum  pen  recorder.  This  same  device  has 
been  used  in  previously  published  work  (Scott, 
1970).  The  device  could  be  fined  with  a conical-nose 
penetrometer  probe  or  a gravity  cover.  Details  of  the 
tests  and  analysis  of  the  data  have  been  presented  by 
Mjgliote  and  Lee  (1971),  along  with  a proposed  pre- 
diction method.  The  use  of  this  prediction  method 
will  be  discussed  further  in  a following  section. 

Large  Seafloor  Penetrometer 

In  an  effort  to  obta  n data  to  verify  the  scaling 
relationships  derived  from  the  model  tests,  full-scale 
tests  were  conducted  with  a large  penetrometer  in 
terrigenous  silt  in  the  Santa  Barbara  Channel.  The 
device  is  shown  in  Figure  16;  it  is  rcady  for  drop  from 
a surface  ship  in  1,200  feet  of  water.  The  penetro- 
meter is  about  25  feet  iong  and  weighs  about  2,600 
pounds;  it  was  constructed  by  filling  a 6-inch-outside- 
dia meter  steel  tube  with  lead. 

Attempts  were  made  to  instrument  its  nose  so 
that  force  could  be  measured  and  rigid-bedy  decelera- 
tion recorded  with  an  on-board  accelerometer.  The 
instrumentation  system  performed  poorly,  yielding 
no  force  data  and  only  one  readable  trace  of  decelera- 
tion for  the  several  tests  conducted.  The  test  data  arc 
summarized  in  Table-  2.  The  observed  performance 
will  be  compared  with  other  data  in  a following 
section. 

Air-Dropped  Projectiles 

Tests  of  models  having  various  shapes  and 
various  projecting  fins  were  conducted  to  provide 
additional  data  on  the  effects  shape  and  projections 
have  on  penetration  resistance.  The  penetrators  were 
assembled  from  interchangeable  parts  designed  so 
that  many  different  combinations  of  shape  and  fin  set 
could  be  obtained.  Enough  parts  were  fabricated  to 
assemble  four  complete  penetrators  at  any  one  time, 
such  as  shown  in  Figure  17. 


Figure  16.  Large  seafloor  penetrometer  being 
launched  from  MV  Gear. 


A site  was  selected  on  the  western  boundary  of 
the  Naval  Shipyard- at  Mare  Island,  California,  near 
the  main  sewage  outfall  for  the  station.  Soil  at  the 
site  is  recently  deposited  San  Francisco  Bay  mud  in  a 
nearly  flat,  straight, shore  area  maintained  naturally  in 
a wet  condition  year-round.  Strength  data  for  soil  at 
the  sewage  outfall  were  available  from  logs  of  soil 
tests  conducted  for  construction  of  the  outfall;  these 
data  arc  shown  in  Figure  18.  The  data  indicate  that 
the  sediment  down  to  about  30  feet  has  a ratio  of 
shear  strength  to  overburden  pressure  of  approxi- 
mately 0.3,  and  that  the  underlying  sediment  is 
substantially  stronger,  possibly  because  of  desiccation 
or  erosion  that  occurred  during  the  last  geological 
recession  of  the  sea.  Based  upon  the  topographical 
uniformity  of  the  site,  these  properties  were  general- 
ized to  the  adjacent  target  area  200  to  300  feet  north 
of  the  outfall. 

A helicopter  was  utilized  as  a work  platform  for 
dropping  and  extracting  the  projectiles.  Each  test 
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Figure  17.  Uninstrumented,  droppable  field  pcnctrators;  parts  arc  interchangeable. 


consisted  of  (1)  hanging  a projectile  50  feet  beneath 
the  aircraft  on  a pullout  line,  (2)  dropping  the  pro- 
jectile from  an  altitude  of  approximately  200  feet  by 
releasing  the  pullout  line,  (3)  hovering  at  an  altitude 
of  20  feet  to  retrieve  the  pullout  line  by  hauling  in  a 
tag  line  attached  to  the  end  of  the  pullout  line,  and 
(4)  extracting  the  projectile  from  the  mud  by 
increasing  aircraft  altitude  with  the  end  of  the  pullout 
line  attached  to  a load-measuring  device  suspended 
from  a frame  in  the  aircraft. 

The  projectiles  were  not  instrumented.  However, 
high-speed  motion  pictures  were  taken  at  1,000 
frames  per  second  for  observing  any  instability  of  the 
projectiles  when  entering  the  mud,  and  elevated-speed 
motion  pictures  wcic  taken  at  48  frames  per  second 
to  provide  data  on  projectile  impact  velocity,  flight 
stability,  and  embedment  depth.  Color  slides  pro- 
vided back-up  data  on  embedment  depth. 


The  film  data  were  reduced  frame-by-frame 
manually  by  scaling  to  obtain  travel  distances  Depth 
of  embedment  was  determined  by  scaling  the  length 
of  pullout  line  above  the  mud  when  the  line  first 
came  taut  during  extraction,  ar.d  subtracting  this 
length  from  the  known  length  of  the  pullout  line. 
Timing  was  determined  from  the  known  film  speed. 

Thirty-six  such  tests  were  conducted.  In  differ- 
ent tests,  the  projectiles  were  assembled  having  three 
different  length-to-diametcr  ratios  (slendernesses), 
tailpieces  with  and  without  stabilizing  fins,  and  nose- 
pieces  with  and  without  streamlining,  with  no  flukes, 
fairing,  and  with  simulated  rectangular  and  circular 
anchor  flukes  (simulated  flukes  did  not  rotate  to  key 
during  extraction).  Combinations  were  varied  to  pro- 
vide a sampling  of  the  variables  which  was  approxi- 
mately balanced  statistically,  with  repetitions  of  some 
configurations. 
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Figure  18.  Shear  strength  from  unconfincd  compression  tests  at  Mare  Island,  CA, 
test  site  (Anon.,  1956). 


The  test  data  are  shown  in  Table  6 along  with 
the  analytical  calculations  discussed  in  a following 
section.  The  tests  are  grouped  according  to  slender- 
ness. 


Fgg  = bearing  pressure  force 


AD 


side  adhesion  force 


fh  = fluid  inertial  drag  force 


DATA  ANALYSIS 
Underlying  Laws 

The  physical  laws  pertaining  to  penetration  resis- 
tance are  related  to  the  forces  shown  in  Figure  19. 
Vertical  force  equilibrium  gives 

Net  Force  = F„  + W’  - F (la) 

F = FBE  + FAD  + fH  <lb> 

where  FD  = externally  applied  driving  force 

W'  = buoyant  weight  of  penetrator  in  soil 
F = total  penetration  resistance  force 


The  buoyant  weight  of  the  penetrator  and  the 
externally  applied  driving  force  contribute  to  penetra- 
tion, whereas  the  point  bearing  pressure  force,  side 
adhesion  force,  and  fluid  inertial  drag  force  resist 
penetration.  The  direction  of  all  vector  quantities  in 
this  figure  is  positive  downward.  These  forces  consti- 
tute the  analytical  framework  used  in  this  study-,  the 
relative  magnitudes  and  underlying  mechanisms  for 
these  forces  were  left  open  to  adjustment  according 
to  the  data.  Any  net  unbalance  in  these  forces  would 
cause  a change  in  the  penetrator  velocity  by  accelera- 
ting the  mass  of  the  penetrator  along  with  the  added 
mass  of  the  soil  moving  in  association  with  the  penc- 
trator.  The  conditions  of  zero  driving  force,  a 
negative  penetration  resistance  force,  and  a positive 
initial  velocity  correspond  to  free-fall  penetration. 
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Table  6.  Summary  of  Air-Dropped  Projectile  Data 


— 

Ratio  of  Calculated  to  Measured 

Test 

Slenderness 

Ratio, 

L/D 

Initial 

Velocity, 

Measured 

Embedment, 

Depth  of  Embedment  for  a 
Side  Minimum  Adhesion 

No. 

voJ 

z 

Factor,  6*.  of- 

(ft/sec) 

(ft) 

- 



' 

0 

0.1 

0.2 

0.3 

0.4 

0.5 

ADS 

9.3 

60* 

14 

1.21 

1.11 

1.04 

AD10* 

9.3 

66 

17 

1.08 

1.01 

0.96 

AD15 

9.3 

60* 

11 

1.55 

1.44 

1.34 

AD16 

9.3 

60* 

16 

1.02 

0.93 

0.87 

AD20 

9.3 

51 

14* 

1.09 

1.01 

0.96 

AD21* 

9.3 

50 

14 

1.07 

0.98 

0.91 

AD22 

9.3 

60* 

24 

0.69 

0.63 

0.59 

AD27 

9.3 

60* 

13 

1.27 

1.16 

1.08 

AD32 

9.3 

60* 

12 

1.43 

1.30 

1.21 

AD34 

9.3 

61 

18‘* 

0.95 

0.86 

0.80 

AD36 

9.3 

65 

U<' 

1.53 

1.42 

1.32 

Average 

9.3 

1.15 

1.06 

0.99 

Selected 

Average* 

9.3 

1.08 

1.00 

0.94 

AD9C 

15.5 

63 

19 

0.96 

0.88 

15.5 

65 

16rf 

1.24 

1.11 

35.5 

63* 

19 

0.92 

0.85 

AD17 

»5;5 

56 

20J 

0.89 

0.81 

AD19 

15.5 

60 

\5-' 

1.12 

1.03 

AD24* 

15.5 

66 

17 

1.09 

0.99 

AD25 

15.5 

67 

17  d 

1.08 

0.99 

AD28 

15.5 

63* 

14 

1.31 

1.20 

AD29 

15.5 

63* 

20 

0.93 

0.85 

AD35 

15.5 

57 

,3./ 

1.35 

1.22 

Average 

15.5 

1.09 

0.99 

■ 

Selected 

Average* 

15.5 

1.03 

0.94 

continued 
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Tabic  6.  Continued 


Ratio  of  Calculated  to  Measured 

Slenderness 

Ratio, 

L/D 

Initial 

Measured 

Depth  of  Embedment  for  a 

Test 

Velocity, 

Embedment, 

Minimum  Side  Adhesion 

No. 

v a 
vo 

(ft/sec) 

z 

(ft) 

Factor,  6*,  of- 

■■■ 

D 

0.1 

mm 

0.3 

0.4 

1 

AD4 

■ 

67* 

23 

0:82 

0.76 

ADU 

65 

1SJ 

0.99 

0.92 

AD13C 

27.8 

72 

21 

1.05 

0.96 

AD18C 

27.8 

60 

17 

1.03 

0.94 

AD23 

27.8 

56 

lO** 

1.64 

1.51 

AD26 

27.8 

66 

7d 

2.58 

2.38 

27.8 

67* 

22 

0.85 

0.79 

AD31 

27.8 

67* 

UJ 

1.67 

1.34 

AD33 

27.8 

67* 

17 

1.13 

1.04 

Average 

27.8 

■ 

'/ 

m 

1.07 

0.99 

Selected 

Averagec 

27.8 

■ 

■ 

■ 

1.04 

0.95 

a Reduced  from  impact  velocity  to  account  for  effects  of  energy  transfer  to 
added  mass  on  impact. 

* Velocity  measurement  not  obtained;  indicated  value  is  average  for  this 
projectile  shape. 

c Selected  averages  are  for  test  numbers  indicated. 
l/  Questionable  embedment  depth  because  of  test  anomalies. 
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fluid  inertial  drag  force,  FH 


bearing  pressure  force,  Fgg 

Figure  19.  Forces  acting  on  a penetrating  object. 


For  ease  in  analysis,  it  was  assumed  that  the 
fluid  inertial  drag  obeyed  a velocity-squared  law 
identical  to  that  for  a Newtonian  fluid  with  a mass 
density  equal  to  the  bulk  mass  density  of  the  soil 
(Robertson  and  Pazwash,  1971).  The  data  obtained  in 
this  study  fell  in  a range  where  the  soil  strength  over- 
shadowed the  drag  to  the  extent  that  it  was 
impossible  to  verify  this  assumption;  however,  satis- 
factory agreement  is  apparent  and  the  assumption  is 
not  disputed. 

A formula  that  is  generally  used  to  predict  pene- 
tration behavior,  which  accounts  for  the  forces 
mentioned  above,  is  the  Poncelct  Equation  (Schmid, 
1969),  wherein  the  total  resistance  to  penetration  is 

F = F*  + F„(v2)  (2a) 

where  F¥  = plastic  resistance  force,  a constant 
not  dependent  on  velocity  or 
displacement  (represents  FBE 
and  Fy^u  in  Equation  lb) 

F„(v2)  = fluid  inertial  drag,  proportional 
to  the  square  of  the  velocity 

The  fluid  inertial  drag  may  be  expressed  in  a manner 
commonly  used  for  fluid  flow, 


F„(v2>  = Y^CoApV2  (2b) 

where  p = fluid  mass  density 

CD  = fluid  drag  coefficient  based  upon 
frontal  area 

Ap  = penetrator  frontal  area 
v = penetrator  velocity 

Equation  2a  was  modified  rather  substantially  by 
replacing  the  constant  F*  term  with  a term  depen- 
dent on  both  depth  and  velocity  in  order  to  account 
for  the  nonconstant  nature  of  the  plastic  resistance 
forces  observed  in  the  model  tests.  The  bearing  pres- 
sure force  depends  on  soil  strength,  depth,  and 
patterns  of  soil  deformation.  The  side  adhesion  force 
depends  on  soil  strength,  reduction  in  soil  strength  by 
remolding,  and  adhesion  or  nonadhesion  (separation) 
between  the  side  of  the  penetrator  and  the  adjacent 
soil.  Further,  the  soil  strength  normally  varies  with 
depth  and  is  known  to  depend  on  strain  rate,  which  is 
related  to  the  velocity  and  the  geometry  of  the  penc- 
trator.  A modified  form  having  features  to  account 
for  these  effects  was  hypothesized,  incorporating 

F*  * Si<FBE  + FAD>  (2c) 
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where  S*  = soil  strain-rate  effect,  the 

ratio  of  soil  strength  at  a 
velocity  to  that  at  zero 
velocity 

FB8  and  fad  = bearing  and  adhesion  forces, 
respectively,  dependent  on 
static  soil  strength,  depth, 
and  geometry 

The  bearing  pressure  and  side  adhesion  forces  are 
given  by  the  conventional  formula  (Tcrzaghi  and 
Peck,  1967,  or  other  standard  soils  text): 

fBE  = SuNcAp  (3) 


in  Tables  3 through  S,  along  with  the  subsequently 
discussed  calculated  parameters. 

The  instantaneous  penetration  data  were  norm- 
alized and  plotted  in  an  effort  to  discover  significant 
trends  of  variation  in  the  dynamic  soil  strength  and 
the  amount  of  side  adhesion  occurring  at  various 
velocities  throughout  the  penetration  events.  A 
measure  of  the  dynamic-to-static-soil-strength* 
(termed  the  strain-rate  effect)  ratio  was  obtained  by 
subtracting  the  computed  fluid  inertial  drag 
(Equation  2b)  from  the  total  penetration  resistance 
(product  of  penetrator  mass  and  measured  decelera- 
tion, augmented  by  penetrator  buoyant  weight),  and 
dividing  this  difference  by  a sum  of  products  of  areas 
and  resistance  factors.  The  resulting  expression  is 


‘AD 


(4) 


where 


Su  = static  undrained  soil  shear  strength 

Nc  = bearing  capacity  factor 

AF  = penetrator  frontal  area 

5 = side  adhesion  factor 

As  = penetrator  side  area 

St  = soil  sensitivity,  ratio  of  undisturbed 
to  remolded  static  shear  strengths 


The  similarity  between  the  dynamic  deformation 
patterns  observed  in  the  high-speed  movies  (Figure 
13)  and  the  deformation  patterns  associated  with 
conventional  static  bearing  capacity  factors  led  to  the 
adoption  of  conventional  static  bearing  capacity 
factors  for  the  Nc  term.  The  strain-rate  effect,  S^,  and 
the  side  adhesion  factor,  6, remained  as  undetermined 
functions  of  penetrator  velocity  and  geometry,  and 
soil  properties. 


where  F*  = F - F„.  This  is  effectively  the  ratio 
of  dynamic  measured  strength  to  static  predicted 
strength. 

Values  for  this  ratio  were  calculated  from  the 
available  data;  a pair  of  values  was  calculated  from 
each  data  set  for  assumed  values  of  side  adhesion 
factor,  8,  of  zero  and  unity.  This  procedure  provided 
a range  of  possible  values  of  the  measured  strain-rate 
effect  ratio  for  each  data  set  corresponding  to  values 
of  the  side  adhesion  factor  ranging  between  the 
assumed  extremes.  The  exact  values  for  the  strain-rate 
effect  and  the  side  adhesion  factor  existing  under  any 
instantaneous  set  of  circumstances  were  hypothesized 
to  depend  on  velocity,  soil  properties,  and  penetrator 
characteristics. 

In  order  to  determine  the  nature  of  the  varia- 
tions in  the  strain-rate  effect  with  velocity  and  other 
conditions,  a normalized  velocity  was  selected  as 


Point-By-Point  Instantaneous  Data 

Several  points  were  selected  at  equal  intervals  in 
time  for  each  instrumented  model  test  in  the  silt- 
cement  test  series  and,  with  the  exception  of  the 
half-section  testa,  for  each  test  in  the  soft  and 
medium  silt  series.  The  instantaneous  deceleration, 
velocity,  and  displacement  for  these  points  arc  given 


where  L is  the  effective  length  of  the  shearing  zone, 
equal  to  the  length  of  the  penetrator  body.  This  term, 
when  multiplied  by  a dynamic  viscosity,  reflects  the 
ratio  of  viscous  forces  to  plastic  forces  in  a flowing 
fluid.  This  term  was  found  to  correlate  well  with 


• "Static"  was  taken  as  denoting  behavior  at  a strain  rate  equal  to  that  employed  in 
strength  testing  of  soft  silt. 
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variations  both  in  the  strain-rate  effect  and  in  the  side 
adhesion  factor.  Instantaneous  data  are  shown  on 
plots  of  computed  ranges  of  S*  tersus  v*  in  Figures 
20  through  22.  In  these  plots;  the  S£  ranges  are 
shown  corresponding  to  values  of  8 between  the 
extremes  zero  and  unity.  It  was  anticipated  that,  as 
velocity  increased  from  zero  to  a high  value,  S£  would 
increase  from  1.0  to  a somewhat  higher  value,  and  8 
would  decrease  within  the  extremes  unity  and  zero. 
Thus,  a function  was  developed  to  represent  the 
strain-rate  effect  at  all  values  of  v*.  Finite  values  of 
S-  were  maintained  at  the  extremes  in  velocity  of 
zero  and  infinity.  The  developed  function  is 


_L  = i i 

Si  <C£v*  + C0),/2  S«* 


(7) 


C0  « 0.04  (9c) 

It  is  shown  below  that  the  revised  curve  also  fits  the 
overall  data  as  well  as  the  earlier  curve. 

Values  of  5 are  indicated  by  the  position  of  the 
strain-rate  effect  curve  on  the  plotted  ranges  in 
Figures  20  through  22.  As  shown,  8 tends  to  be 
nearly  unity  at  lower  values  of  the  abscissa  for 
pointed  penetrators  in  silt,  and  approximately  0.8  to 
0.9  for  blunt  penetrators  under  the  same  conditions. 
At  higher  velocities,  8 drops  to  vaguely  defined  values 
between  0.0  and  0.5.  This  observed  behavior  agrees 
well  with  the  conclusions  drawn  by  True  (1974)  for 
the  overall  data; 

Overall  Data 


where  S£*  = maximum  S£  at  high  velocities, 

referenced  to  the  normalized  velocity 
of  2 x 10"4  m2/kgf-sec  employed  in 
the  static  strength  tests  of  soft  silt* 

C"  = constant,  having  units  of  stress  times 
time 

C„  = constant,  dimensionless 

In  an  earlier  publication  (True,  1974),  the  three 
arbitrary  constants  in  this  function  were  evaluated  for 
overall  conditions  (impact  velocity,  ultimate  embed- 
ment depth)  as 


S£*  = 3.7 

(8a) 

_ „„„  kgf-scc 

Ci  = 200  — 

c 2 

(8b) 

m 

C = 0.08 

(8c) 

however,  trial  fits  with  the  instantaneous  data  showed 
that  somewhat  different  values  yielded  a substantially 
better  fit  for  the  instantaneous  data;  they  are 


S-*  = 5.0 

(9a) 

« kgf-scc 
= 100  -2 

(9b) 

m2 

The  overall  data  and  the  relationships  derived 
therefrom  have  been  presented  previously  by  True 
(1974).  The  overall  data  refers  to  measurements  of 
overall  behavior,  including  impact  velocity  and  ulti- 
mate embedment  depth,  and  computed  parameters 
based  on  these  measurements.  The  information 
presented  here  deviates  somewhat  from  that  previous 
presentation,  and  it  relates  the  ‘‘overall”  penetration 
performance  to  the  instantaneous  behavior.  The  over- 
all data  from  the  model  tests  are  presented  in  Table  2 
along  with  the  major  computations  in  tabular  form 
leading  to  the  graphical  derivation  of  the  penetration 
resistance  relationships.  The  half-section  tests  have 
been  omitted  as  the  results  from  these  tests  are 
distorted  by  the  effects  of  friction  between  the  pene- 
trator  and  the  transparent  plastic  face  of  the  test 
tank. 

The  final  columns  in  Table  2 show  ranges  of 
strain-rate  effect  calculated  from  the  presented  data 
for  the  limiting  values  of  adhesion  reduction  factor, 
zero,  and  unity.  These  ranges  are  shown  plotted 
against  normalized  velocity  in  Figure  23.  The  line 
represented  by  Equation  7 is  also  plotted,  with  the 
constants  as  given  by  Equations  9a  through  9c.  This 
line  follows  the  lower  ends  of  the  plotted  ranges  of  S£ 
(corresponding  to  5 = 1)  for  low  velocities,  transiting 
to  the  upper  ends  of  the  range  (corresponding  to  8 = 
0)  for  higher  velocities.  The  relationship  from  the 
earlier  reported  analysis  (True,  1974)  of  the  overall 
data,  as  given  by  Equations  7 through  8c,  is  shown  as 
a dotted  curve  for  comparison.  It  is  apparent  that  the 


* Static  strengths  measured  at  other  values  of  normalized  velocity  require  normalization 
to  this  value. 
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revned  (solid)  curve,  as  derived  herein  from  the 
analysis  of  the  data  on  an  instantaneous  basis,  fits  the 
overall  data  as  well  as  the  previously  derived  curve. 

Values  of  8 were  back-calculated  from  the  inter- 
sections of  die  fitted  (solid)  line  with  the  plotted 
ranges;  these  arc  shown  plotted  against  the  norma- 
lized velocity  in  Figure  24.  From  this  plot  it  is 
apparent  that  8 assumes  values  near  0.8  at  low 
velocities,  transiting  to  low  values  at  high  velocities, 
with  transition  zones  and  high-velocity  values  of  5 
varying  in  position  for  different  nose  shapes. 

The  values  of  6 at  high  velocities  and  the 
positions  of  the  transition  zones  should  depend  on 
the  length-to-breadth  (slenderness)  ratios  of  the  penc- 
trators.  As  most  of  the  penetrator  travel  in  dynamic 
penetration  occurs  at  velocities  above  the  transition 
velocity,  the  value  of  8 at  high  velocities  is  of  primary 
concern.  It  is  reasonable  that  this  value  of  8,  termed 
8*,  for  more  slender  penetrators  should  be  greater 
than  or  equal  to  die  predicted  values,  and  for  the 
stubbier  penetrators  it  should  be  less  than  or  equal  to 
the  predicted  values  based  upon  the  present  model 
tests. 

The  field  tests  of  the  penetrators  of  varying 
slenderness  conducted  in  San  Francisco  Bay  mud  at 
Mare  Island  provide  data  for  determining  the  effect  of 
slenderness  on  the  asymptotic  value  of  the  side 
adhesion  factor,  5*,  at  high  velocities.  The  measured 
depths  of  embedment  ranged  from  10  to  25  feet.  The 
calculated  embedment  depths  were  obtained  by 
substituting  trial  values  of  5*  together  with  esti- 
mated values  of  other  parameters  into  an  incre- 
mental calculation  procedure  (described  in  a succeed- 
ing section).  The  values  of  8*  were  adjusted  until  the 
calculated  embedment  depths  were  in  best  agreement 
with  the  measurements,  as  shown  in  Table  6.  It  was 
found  that,  by  making  a suitable  selection  of  the  5 * 
value,  discrepancies  between  calculated  and  measured 
values  of  embedment  depth  could  be  reduced  to  a 
level  consistent  with  natural  variations  in  soil 
strength.  The  resulting  values  of  6*  are  summarized  in 
Figure  25  and  Table  7.  As  the  embedment  depths 
computed  by  the  incremental  procedure  vary  only 
about  10%  for  variations  in  6*  of  over  20%,  the  selec- 
tion of  approximate  values  of  5*,  using  Table  8 as  a 
guide,  should  yield  sufficient  accuracy  for  most  pur- 
poses. 


Table  7.  Minimum  Values  of  High-Velocity 
Side  Adhesion  Factor  Derived 
From  Field  Test  Data 
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The  field  tests  of  the  large  seafloor  penetrometer 
provided  data  for  checking  the  above-derived 
relationships  for  a large,  slender  penetrator.  The  test 
data  are  summarized  in  Table  2.  The  back-calculated 
values  of  6 are  greater  than  unity;  possible  causes  for 
this  discrepancy  include  the  effect  of  an  unac- 
counted-for line  drag  force,  erroneously  low  values  of 
assumed  soil  strength  at  sub-sampled  depths,  erro- 
neously low  values  of  recommended  soil  strain-rate 
effect  at  the  appropriate  value  of  velocity  term,  or  a 
velocity  term  which  does  not  properly  represent  wide 
ranges  in  penetrator  size  and,  especially,  penetrator 
length  as  related  to  the  length  of  shearing  path.  The 
data  are  insufficient  to  resolve  this  question. 

DISCUSSION 

Components  of  Penetration  Resistance 

The  resistance  imparted  by  a soft,  saturated  soil 
to  the  deep,  vertical  penetration  of  an  axisymetric 
cylindrical  penetrator  can  be  resolved  into  com- 
ponents in  order  to  enable  generalization  to  soils  of 
various  strengths  penetrated  at  various  velocities.  To 
this  end,  the  observations  made  in  the  present  study 
were  considered  in  light  of  previously  established 
theory  and  empirical  information. 

Bearing  Pressure 

Observations  of  the  deformation  patterns  of  soil 
around  an  advancing  half-section  penetrator  showed 
no  separation  or  other  visible  departure  from  patterns 
observed  after  penetration.  In  all  cases,  these  patterns 
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appeared  to  duplicate  patterns  associated  with  the 
development  of  conventional  static-bearing  capacity 
factors  based  on  plasticity  theory.  This  led  to  the  con- 
clusion that  the  conventional  static  bearing  pressure 
theory  is  applicable,  in  a general  sense,  to  the  phe- 
nomenon of  dynamic  penetration,  subject  to  modifi- 
cation only  to  the  extent  necessary  to  account  for  the 
inertial  forces  in  the  system  and  the  effects  of  high 
rates  of  shear  strain  on  the  shearing  resistance  in  the 
soil.  Inertial  forces  are  associated  with  the  rigid-body 
accelerations  of  the  pcnetrator  and  its  effective  added 
mass  and  the  fluid-type  inertial  drag  in  the  soil.  Vis- 
cous drag  forces  are  lumped  with  the  soil  shearing 
resistance  in  this  analysis. 

Side  Adhesion  Factor 

The  adhesive  shear  stress  along  the  sidewall  of 
the  pcnetrator  normally  is  considered  separately  from 
the  bearing  pressure  in  conventional  pile  capacity 
analysis.  The  adhesive  shear  stress  immediately  after 
driving  takes  on  a value  equal  to  the  product  of  the 
remolded  strength  of  the  soil  multiplied  by  a side 
adhesion  factor.  The  side  adhesion  factor,  normally 
unity  for  soft  soils,  was  suspected  of  being  subject  to 
reduction  by  the  tendency  of  the  soil  to  separate 
from  the  penetrator  sidewall.  This  tendency  was 
hypothesized  to  result  from  the  outward  momentum 
of  the  soil  diverted  around  the  nose  of  the  advancing 
penetrator.  Such  conditions  would  tend  to  reduce 
effective  soil  stresses  and  increase  pore  pressures, 
thereby  reducing  the  net  soil  strength  adjacent  to  the 
penetrator.  The  increase  in  pore  pressures  would  have 
been  facilitated  by  previously  supernatant  water 
tending  to  be  carried  along  with  the  penetrator  as  it 
entered  the  soil.  Thus,  substantial  reductions  in  the 
side  adhesion  factor  were  anticipated  at  high  veloci- 
ties, even  in  very  soft  soils. 

The  values  given  ir.  Table  7 are  recommended  as 
a guide  in  selecting  a value  of  6 * for  a penetrator  of 
concern.  If  more  precision  is  required,  model  tests 
can  be  conducted  to  obtain  measurements  of  the 
adhesion  factor  for  any  shape  of  interest. 

Strain-Rate  Effect 

Sufficient  data  were  obtained  on  the  soil 
materials  used  In  this  work  to  determine  the  effect  of 
strain  rate  on  the  shearing  resistance  of  the  silt  test 


soil.  The  vane  shear  test  data  were  obtained  at  rates 
of  ?3, 166,  and  830  degrees  per  minute  on  the  softer 
model  test  materials,  and  unconfined  and  confined 
compression  test  data  were  obtained  on  the  stronger, 
cemented  materials.  Furthermore,  the  penetration 
tests,  which  were  conducted  at  various  rates,  provided 
substantial  additional  information  on  strain-rate 
effect. 

The  phenomenon  of  penetration  offers  some 
distinct  advantages  as  a basis  for  determining  the 
effect  of  strain  rate  on  shearing  resistance  in  a soil.  It 
can  be  conducted  at  high  rates  with  inertial  effects 
limited  to  a relatively  small  fluid  inertial  drag  force. 
The  tendency  toward  separation  at  the  sidewall  can 
be  controlled  by  properly  adjusting  die  nose  shape 
and  breadth  of  the  afterbody  and/or  it  can  be  moni- 
tored with  appropriate  instrumentation.  Finally,  as 
no  significant  variation  has  been  found  in  the  shape 
of  deformation  patterns  around  advancing  penc- 
trators  at  varying  velocities,  the  use  of  a constant 
bearing  pressure  coefficient  should  yield  an  accurate 
picture  of  variations  in  soil  shearing  resistance  with 
strain  rate  when  applied  to  penetration  data  obtained 
at  different  velocities. 

The  problem  remains  of  relating  the  velocity 
difference  between  two  blocks  of  soil  or  between  a 
foreign  body  penetrating  soil  to  the  shear  strain  rate. 
If  a constant  deformation  pattern  is  maintained,  it  is 
reasonable  to  assume  that  the  shear  strain  rate  is 
nearly  proportional  to  the  velocity.  Analogous  to  the 
behavior  of  a fluid  boundary  layer,  it  is  probable  that 
the  thickness  of  a shearing  zone  is  proportional  to  its 
length.  As  the  magnitude  of  shear  strain  for  a given 
overall  displacement  varies  inversely  with  this 
thickness,  it  is  reasonable  to  assume  further  that  the 
shear  strain  rate  is  inversely  proportional  to  the 
length  of  the  shearing  zone.  This  length  corresponds 
to  the  circumference  of  the  vane  in  a vane  shear  test 
and  roughly  to  the  length  of  the  penetrator  in  the 
case  of  penetration.  These  two  assumptions  appear 
reasonable  and  were  adopted  for  the  present  analysis. 
The  resulting  derived  relationships  are  recommended 
for  engineering  applications.  For  purposes  of  pre- 
dicting penetrator  performance,  the  relationship 
derived  above  for  the  strain-rate  effect  (Equations  7, 
9a,  9b,  and  9c)  should  be  used  to  compute  a dynamic 
soil  strength  in  soils  composed  of  strong  grains 
relative  to  the  bulk  soil  strength.  This  dynamic  soil 
strength,  i.e.,  the  product  of  strain-rate  effect  and 
undrained  strength,  is  to  be  utilized  in  calculating  the 
soil  resistance  force. 
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Soil  Shear  Strength 


Other  Parameters 


The  undrained  shear  strength  of  a target  material 
can  be  determined  by  conventional  sampling  and 
testing  methods  normally  employed  in-'he  course  of 
any  sound  foundation  engineering  effort.  However,  it 
often  is  use  ful  to  employ  a genet.;/  ruk  of  ‘.Sumb  as 
an  expedient  in  the  preliminary  stages  of  an  effort. 
For  clays,  soft  silts,  and  some  finegrained  seafloor 
sediments  comprising  a large  portion  of  the  area  of 
the  seafloor,  such  a rule  of  thumb  can  be  surprisingly 
accurate.  In  general,  the  undrained  shear  strength  in 
the  upper  few  tens  of  feet  in  such  sediments  increases 
approximately  linearly  with  depth  from  essentially 
zero  at  the  sediment  surface,  provided  the  sedimenta- 
tion process  is  relatively  constant  with  no  removal  of 
upper  layers  by  erosion,  slumping,  or  other  means. 
The  constant  of  proportionality’,  termed  C,,  has  units 
of  stress  divided  by  length  (depth).  Values  of  Cj 
depend  upon  the  soil  composition  and  density,  falling 
mainly  in  the  range  between  80  and  240  (kgf/m2)/m. 
Without  more  specific  information,  a value  may  be 
assumed  for  Cz  for  a generalized  silt-sized  or  finer 
deep  seafloor  sediment  as 


Cz 


160  k£f/m.2 
m 


(10) 


A constant  strength  term  can  be  added  to  the  product 
of  Cz  with  depth  if  an  elevated  overall  strength  is 
anticipated. 

Similarly,  the  ratio  of  undisturbed  to  remolded 
strength  for  fine-grained  seafloor  sediments,  termed 
the  sensitivity,  St,  varies  from  approximately  2 for 
medium  silts  to  approximately  5 for  clays;  a value  can 
be  assumed  for  St  for  a generalized  silt-sized  or  finer 
deep  seafloor  sediment  as 

St  = 3 (11) 


Based  upon  the  data  obtained  in  this  study,  an 
added  mass  of  twice  the  mass  of  soil  displaced  by  the 
penetrator  accounts  for  the  effects  of  instantaneous 
deceleration  on  observed  behavior;  this  value  of 
added  mass  is  recommended  for  use  in  dynamic  pene- 
tration predictions. 

This  study  showed  no  reason  to  doubt  the  use  of 
standard  methods  for  evaluating  the  fluid  inertial  drag 
component  of  penetration  resistance.  However,  as  the 
study  was  conducted  under  conditions  where  the 
magnitude  of  the  drag  component  was  small  relative 
to  the  total  resistance,  this  approach  cannot  be  valida- 
ted by  the  data  obtained.  The  use  of  standard 
methods  for  the  evaluation  of  fluid  inertial  drag  on 
dynamically  penetrating  bodies  is  recommended 
pending  possible  future  developments  to  the  con- 
trary. 

Equation  of  Motion 

In  mathematical  terms,  a penetration  resistance 
force  composed  of  the  components  outlined  above  is 
expressed  by  Equations  2a  through  2c,  3,  and  4, 
which  can  be  combined  to  give 


F 


s«s; 


ncaf 


+ \p  CDApV2  (12) 

From  Equation  la  and  Newton’s  second  law  of 
motion,  the  equation  of  motion  is 

M*  4~  = w*  + FD  - F (13) 


It  is  emphasized  that  these  approximations  are  for 
preliminary  efforts  only  and  should  be  replaced  by 
more  definite  representations  of  soil  properties  when 
making  more  precise  engineering  predictions. 


where 


M*  = penetrator  mass  plus  added  mass 
t = time 

d = differential  operator 
W*  = penetrator  buoyant  weight 
Fd  = penetration  driving  force 
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CALCULATION  PROCEDURE 
Development  of  Procedure 

The  equation  of  motion  for  a penetrator  formu- 
lated to  account  for  the  effects  of  strain  rate,  side 
adhesion,  and  other  factors  discussed  in  the  preceding 
section  is  not  suitable  for  solution  in  closed  form. 
However,  it  can  be  solved  readily  by  incremental 
techniques.  The  purpose  of  this  section  is  to  present 
such  a technique  by  which  penetration  behavior  can 
be  estimated  relatively  precisely  with  reasonable 
hand-calculation  effort  or  with  only  seconds  of 
calculation  time  on  a desk-top  computer. 

The  equation  of  motion  for  a penetrator  is  given 
by  Equation  13,  with  penetration  resistance  force 
given  by  Equation  12.  A different  form  of  the  equa- 
tion of  motion,  with  time  eliminated  to  facilitate 
hand  computations,  is 

M*v^  = W'  + F„  - F(v, z)  (14) 

The  incremental  form  to  be  used  for  hand  computa- 
tions is 

lW,+  FDj-F(vi,2iH 

vi+l  - vi-l  + — m^: <*i+l  - <15> 

In  this  form,  all  functions  are  known  except  v;+ 1 ; 
(Zj+I  - Zj.j)  is  specified  to  be  one-tenth  or  less  of  a 
guessed  embedment  depth.  When  beginning,  it  is 
necessary  to  guess  Vj ; this  is  done  most  directly  by 
computing  v2  for  vt  = v0  and  then  starting  over  again 
using 


It  is  not  necessary  to  guess  Vj  more  precisely  than 
this  to  yield  a near-perfect  final  value,  zn,  as  vn 
approaches  zero.  A flow  diagram  for  this  technique 
appears  in  Figure  26. 

Comparison  with  Other  Prediction  Methods 

The  predictive  accuracy  of  the  relationships 
derived  herein  was  compared  to  that  of  previously 
available  methods.  A broad  range  of  penetrator  sizes 


and  associated  penetration  distances  was  used. 
Included  were  a test  of  the  small,  blunt  steel  model  in 
a soft  silt  target,  an  air-drop  test  of  the  medium- 
slenderness  model  in  the  natural  shoreline  deposit  of 
soft  silty  clay  at  Mare  Island,  and  a test  of  the  large 
seafloor  penetrometer  in  a soft  terrigenous  clayey  silt 
laying  beneath  1,200  feet  of  water  in  the  Santa 
Barbara  Channel.  The  depths  of  penetration  for  these 
tests  were  measured  at  2.3  feet,  17  feet,  and  23  feet, 
respectively.  Targets  were  limited  to  soft  cohesive 
soils.  “Predicted”  depths  of  penetration  were  calcu- 
lated by  several  methods  for  the  documented 
conditions  occurring  during  the  selected  tests;  these 
“predictions”  were  compared  with  the  test 
measurements  of  penetration  depth.  The  resuits  are 
plotted  in  Figure  27. 

From  this  plot  it  is  apparent  that  the  present 
work,  the  method  of  Migliore  and  Lee  (1971),  and  a 
method  involving  a Poncelet-form  equation  with  a 
static  pile  bearing  capacity  used  as  the  constant  com- 
ponent of  penetration  resistance  give  reasonable 
results;  the  largest  discrepancies  between  observed 
and  predicted  depths  were  for  the  large,  seafloor 
penetrometer  for  which  the  present  work  was  the 
most  accurate.  The  methods  of  Young  (1969)  and 
Petry  (as  given  by  Young,  1969,  and  Christians  and 
Meisburger,  1967)  were  developed  primarily  for 
harder  materials  and  yielded  poor  results  when 
extrapolated  to  low  velocities  and  soft  materials  with- 
out empirical  modifications. 

The  method  developed  herein  is  recommended 
as  it  is  more  fundamentally  sound  and  hence  is 
expected  to  be  more  accurate  than  even  the  best  of 
the  other  methods  (Migliore  and  Lee),  particularly  for 
cases  not  closely  duplicated  by  the  experimental  con- 
ditions involved  in  its  derivation.  In  particular,  the 
side  resistance  is  weighted  too  heavily  and  the  frontal 
resistance  is  weighted  too  lightly  in  the  method  of 
Migliore  and  Lee,  and  the  effects  of  remolding  and 
strain  rate  are  not  considered.  Thus,  whereas  the 
method  of  Migliore  and  I.ce  works  well  for  pcnc- 
trators  having  slenderness  ratios  on  the  order  of  5 to 
10  penetrating  soil  having  a sensitiyity  of 
approximately  2 (the  experimental  conditions  upon 
which  it  was  based),  it  is  less  accurate  for  more 
slender  penetrators  (sec  data  for  LSP3  in  Figure  27), 
and  also  is  expected  to  be  less  accurate  for  pene- 
trators of  normal  slenderness  in  target  materials  of 
greater  sensitivity. 
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Example  Problem 

As  an  example,  calculations  are  presented  here 
for  a hypothetical  penctrator,  weighing  300  pounds 
(mostly  lead),  and  configured  to  attain  a relatively 
large  penetration  in  a soft,  seafloor  sediment.  Initial 
setup  calculations  are  completed  to  provide  coeffi- 
cients necessary  to  perform  the  iterative  calculations 
in  a tabular  format.  All  of  the  calculations  generally 
follow  the  flow  diagram  in  Figure  26,  with  slight 
modification  in  ordering  to  facilitate  hand  computa- 
tion. 

The  setup  calculations  are  carried  out  in 
accordance  with  Figure  26  as  follows: 


Az  * 5 ft  (selected  st  about  one-tenth  of 
the  anticipated  ultimate  depth  of 
embedment) 

2.  The  constants  as  defined  in  Figure  26  become 

_ 0.3  \ 2x3x32.2 

Cj  fc  10(9x0.1  ♦ — -x 71 x 5.0 

- 66.3fc/*ec2 


„ 2 x 5 x 32.2  1 96.6 

Cj  ■ s; x — x trr  * 0.7  x 0.1  * 0.007 


307 


32.2 


„ 2x3x237x324 

C3 307 


213ft2/«ec2 


1.  The  parameters  are  established  as  follows: 


C4 


20 

10x3 


0.4  sec 


a.  Material  properties  are  evaluated 

pg  = 96.6  lb/ft3  (typical  for  fine-grained 
seafloor  sediment) 


C,  - 0.04 

3.  The  initial  values  arc 


St  = 3 (Equation  11) 

Cj  = Su/z  = 10  psf/ft  (Equation  10) 

C-  = 20  psf-scc  (Equation  9b) 

C0  = 0.04  (Equation  9c) 

S-*  = 5.0  (Equation  9a) 

b.  Penetrator  characteristics  are 
vr  - 257  lb 

Af  =0.1  ft2 
As  = 7 ft2 

D = 0.33  ft  (penetrator  diameter) 

L = 5 ft 

c.  Penetration  and  flow  parameters  are 

Nc  = 9 (conventional  bearing  capacity  factor) 
CD  = 0.7  (conventional  fluid  drag  coefficient) 
5*  = 0.3  (interpolation  from  Table  8) 

d.  The  initial  velocity  is  stated  as 
v0  = 60  ft/sec 

e.  A calculation  increment  is  selected  as 


Vj  » 60 [0.0674 (60)2  - 21 3J  » 59.2ft/sec 
ij  ■ 5 ft 

*i  -lUr'nb)' 

The  hand-calculated  solution  is  shown  in  tabular 
form  in  Table  8.  The  predicted  depth  of  embedment, 
38  feet,  is  high  but  credible  considering  the  high 
velocity  and  the  low  strength  of  the  deep  seafloor 
sediment. 

The  setup  calculations  described  above  required 
approximately  10  minutes  to  perform;  a few  minutes 
were  required  to  set  up  the  table,  and  the  iterative 
calculations  were  completed  in  another  15  minutes  or 
so.  Thus,  predicting  embedment  depth  by  this  hand- 
calculation  procedure  on  an  occasional  basis  might 
require  a minimum  of  a half  hour  per  test.  By 
utilizing  a desk-top  computer  programmed  to  com- 
plete the  iterative  portion  of  the  computations,  the 
time  required  to  complete  all  computations,  including 
setup,  programming,  and  data  input,  was  reduced  to 
about  10  minutes  per  test  for  the  36  field  tests  of  the 
penetrators  of  varying  slenderness  discussed  in  the 
preceding  section. 
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Table  8.  Example  Hand  Solution  to  Penetrator  Equation  of  Motion 
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FINDINGS  AND  CONCLUSIONS 

1.  The  hypothesized  dependency  of  penetration 
resistance  on  bearing  pressure,  side  adhesion,  and 
fluid  inertial  drag  provides  a satisfactory  framework 
for  predicting  penetration  behavior  in  seafloor 
materials. 

2.  The  effects  of  velocity  and  depth  on  the  bearing 
pressure  and  side  adhesion  components  of  penetra- 
tion resistance  are  significant.  For  deep  penetration, 
they  arc  satisfactorily  accounted  for  by  considering 
their  influence  on  the  dynamic  soil  shear  strength  and 
the  side  adhesion.  For  shallow  penetration,  a reduced 
bearing  capacity  factor  at  shallow  depths  becomes 
significant.  These  effects  arc  given  by  Equations  2a, 
2c,  3,  and  4. 

3.  Model  and  field  test  data  have  shown  that  the 
ratio  of  dynamic  to  static  soil  strength,  termed  the 
strain-rate  effect,  can  be  expressed  as  in  Equation  7 
with  empirical  constants  assuming  values  approxi- 
mately as  given  in  Equations  9a  through  9c.  Simul- 
taneously, data  have  shown  that  the  ratio  of  actually 
developed  to  potentially  available  side  adhesion  resis- 
tance, termed  the  side  adhesion  factor,  varies  from 
near  unity  at  low  velocities  to  some  fraction  at  higher 
velocities,  that  the  high-velocity  value  can  be  used  as 
a constant  for  practical  calculations  of  dynamic  pene- 
tration behavior,  and  that  this  value  varies  approxi- 
mately as  given  in  Table  7 for  blunt-nosed  cylindrical 
pcnctrators  of  varying  slenderness. 

4.  For  the  situations  modeled,  the  strength  of  the 
target  was  sufficiently  high  so  that  the  fluid  inertial 
drag  component  of  penetration  resistance  was  small. 
Hence,  the  hypothesized  velocity-squared 
dependency  is  untested.  It  is,  however,  in  agreement 
with  the  theory  that  the  fluid  inertial  drag  arises  from 
the  change  in  momentum  of  the  soil  as  it  is  deflected 
around  the  advancing  penetrator.  Drag  coefficient 
values,  also  untested  at  this  point,  arc  tentatively 
accepted  as  identical  to  those  for  turbulent  flow  of 
fluids  around  simple  objects  corresponding  to  the 
penetrator  shapes.  A value  of  0.7  is  approximately 
correct  for  cylindrical  penetrators  such  as  those 
tested. 

5.  The  effect  of  acceleration  on  penetration  resis- 
tance was  shown  by  data  on  models  of  different 


masses  to  be  accountable  by  considering  an  added 
mass  along  with  the  penetrator  mass  in  the  equation 
of  motion.  This  added  mass  is  approximately  equal  to 
the  mass  of  twice  the  volume  of  soil  displaced  by  the 
penetrator. 

6.  The  instantaneous  penetration  resistance  force  on 
a penetrator  can  be  estimated  by  using  Equation  12. 
Strain-rate  effect  is  given  by  Equation  7,  with  con- 
stants as  given  by  Equations  9a  through  9c.  The  side 
adhesion  factor  can  be  obtained  by  interpolating  or 
extrapolating  from  the  values  given  in  Table  8;  only 
an  approximate  value  is  required,  as  predictions  are 
not  nearly  as  sensitive  to  side  adhesion  factor  as  to 
strain-rate  effect.  Values  for  soil  strength  and  sensitiv- 
ity can  be  obtained  by  conventional  testing  methods 
or,  for  purposes  of  preliminary  site  evaluation,  can  be 
assumed  to  be  as  given  by  Equations  10  and  11.  A 
bearing  capacity  factor  of  9 should  be  used  for  blunt- 
nosc  cylinders  penetrating  deeply;  for  shallow  pene- 
tration, a reduced  value  should  be  used  in  accordance 
with  conventional  foundation  practice.  A drag  coeffi- 
cient of  0.7  is  satisfactory  for  blunt-nose  cylinders; 
for  other  shapes,  conventional  values  for  turbulent 
fluid  flow  should  be  used. 

7.  The  total  instantaneous  body  force  on  the  penc- 
trator  can  be  used  to  compute  penetrator  motion  and 
ultimate  depth  of  embedment  in  accordance  with 
Newton’s  second  law  of  motion.  A combined  form  of 
the  equation  of  penetrator  motion  suitable  for 
incremental  computation  by  hand  or  on  a desk-top 
computer  is  given  by  Equation  14.  This  form  is  easily 
programmed  for  forward  difference  or  two-sided 
difference  computation. 

8.  The  two-sided  finite  difference  procedure 
described  in  this  report  can  be  used  to  compute  the 
penetration  behavior  of  an  object  into  a sediment  of 
virtually  any  strength-depth  profile.  Computation 
times  are  about  a half  hour  by  hand  or  about  10 
minutes  on  a programmed  desk-top  computer. 

RECOMMENDATION 

A comparison  made  between  several  methods 
for  predicting  embedment  depth  has  shown  that  the 
method  derived  herein  performs  better  than  the 
others  considered.  It  is  recommended  that  this 
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method  be  used  in  the  engineering  of  penetrating 
elements  for  Navy  seafloor  installations. 
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LIST  OF  SYMBOLS 


Instantaneous  deceleration  (L/T2) 

Effective  average  deceleration,  equal  to 
v0z/2Zg(L/T2) 

Penetrator  frontal  area  (L2) 

Pcnetrator  side  area  (L2) 

Fluid  drag  coefficient  based  upon  frontal 
area  (-) 

Normalized  velocity  coefficient  in  strain 
rate  effect  relationship  (FT/L2) 

Constant  in  strain  rate  effect 
relationship  (-) 

Rate  of  increase  in  soil  strength  per  unit 
depth  (F/L3) 

Constants  defined  in  flow  diagram  for 
incremental  calculation  procedure 
(various  units) 

Penetrator  diameter  (L) 

Total  penetration  resistance  force  (F) 

Plastic  resistance  force  of  soil  (F) 

Side  adhesion  force  (F) 

Bearing  pressure  force  (F) 

Externally  applied  driving  force  (F) 

Fluid  inertial  drag  force  (F) 

Combined  net  normalized  resistance  term 
(L2/T2) 

Normalized  soil  resistance  term  (L2/T2) 

Normalized  drag  term  (L2/T2) 

Acceleration  of  gravity  (L/T2) 

Increment  number  in  incremental 
calculation  procedure  (-) 

Penetrator  length  (L) 

Penetrator  mass  (FT2/L) 

Penetrator  effective  mass,  equal  to  mass 
plus  added  mass  (FT2/L) 


Nc  Bearing  capacity  factor  (-) 

n Total  number  cf  increments  in  incremental 

calculation  procedure  (-) 

S;  Soil  strain  rate  effect  (-) 

Sj*  Maximum  value  of  Se  at  high  velocities  (-) 

St  Soil  sensitivity  (-) 

Sa  Static  undrained  soil  shear  strength  (F/L2) 

t Time  (T) 

v Instantaneous  penetration  velocity  (L/T) 

v Penetrator  effective  velocity  for  overall 

event  - equal  to  2v0/3  in  dynamic  tests; 
equal  to  v in  static  tests  (L/T) 

v*  Normalized  velocity  (L2/FT) 

vc  Initial  penetration  entry  velocity  (L/T)  - 

may  be  less  than  incident  velocity  because 
of  loss  of  energy  to  shock  waves  and  to 
acceleration  of  added  mass  during  impact 

W Weight  of  penetrator  ( F) 

W Buoyant  weight  of  penetrator  in  soil  (F) 

w Normalized  weight  term  (L2/T2) 

Z Ultimate  depth  of  embedment  (L) 

z I nstantancous  depth  in  soil  ( L) 

7 Soil  weight  density  (F/L3 ) 

A Finite  incremental  operator  (-) 

6 Side  adhesion  factor  (-) 

6 * Minimum  value  of  6 at  high  velocities  (-) 

p Fluid  (soil)  mass  density  (FT2/L4) 
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